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ON THE CANONICAL FORM OF THE SINGULAR / 


MATRIX PENCIL 
By A. C. AITKEN (Edinburgh) 
[Received 31 July 1933] 


THE object of this note is to derive Kronecker’s canonical form of 
a singular pencil of matrices by elementary rational methods. The 
feature of the method is that the requisite transformations of the 
pencil are studied by way of the reciprocal transformations of certain 
row and column vectors intimately connected with the pencil. 


1. The singular pencil and its minimal vectors 


A pencil of matrices, such as 


A Jel d=— oh) a) 


is singular, in that its determinant vanishes identically in A, ». In 
general a matrix pencil A = AA+B is singular when its determinant 
AA+pB| is identically zero, or when the matrix AA+B is rect- 
angular, not square. 

With respect to example (1) it may be observed that 


fl, —1; “|= 00. 0); (2) 


A p 


also that x 72 = lof (3) 


Indeed the quality of being singular consists in a dependence 
existing among the rows or the columns of the pencil, expressible by 
either or both of the statements 


6A => 0, Ad = 0, (4) 


where @ and ¢ are row and column vectors whose elements are poly- 
nomials in A and p. Since in any element of the vectors 6A and Ad 
all terms of each degree must vanish separately, it follows that there 
exist vectors 6 and ¢ satisfying (4) and such that 


8 = [8,, 0,...,8,], 
d = {1 b+» by} (5) 


where the elements 6; and ¢; are polynomials homogeneous in A, p, 
3695.4 R 








242 A. C. AITKEN 
and each of a certain degree, the same for all the @,, 
for the ¢;. 

Further, as the simple example (1), with (2) and (3), also shows, 


and similarly 


there is a minimal degree for such polynomials, in the sense, for 


A fog 

A pL 
reduces to zero the vector {u,—A}, which has elements of the first 
degree, but not (identically) any vector of the form {1,a}, having 
elements of zero degree. 


instance, that 


We shall use the index m for the minimal degree in vectors 6, and 
m’ for that in vectors ¢; and it will be said that the rows of A have 
dependence of minimal order m in A, », while the columns have de- 
pendence of minimal order m’. 

These minimal orders have invariant properties (3); they are 
invariant under equivalent non-singular transformations of A, and 
also under non-singular linear transformations of A, » to p, o. The 
proof of these statements is immediate; in either case we have only 
to suppose the contrary and then to apply the reciprocal (non- 
singular) transformation, which would yield an identity of less than 
minimal order, involving an absurdity. 

2. Transformation of the minimal vector 

The vector 6 = [6,, 0,,...,4,| (6) 
will now be studied, and transformed. The elements @; are poly- 
nomials of degree m; hence if there are more than m-+-1 non-zero 
elements 6; they cannot be linearly independent. They must be 
expressible linearly in terms of a basic set of m-+1 polynomials, 
perhaps even fewer. We may take as basic set 

pm, —Au™-!, Ayw™-?,..., (—) mar”, (7) 


and we shall prove that none of these can be absent. 


For suppose the contrary. Suppose, for example, that the 6; are 
linearly expressible in terms of m polynomials only, A**4”-*-! being 


missing; so that 
A = [we = Au” “a ye-2 (- . a a k 0. (— ye 2yk i ay-S-8 $ 
| # 


Let a typical column of HA be the column vector 


Y = (QA), p, a,A+bop 








THE SINGULAR MATRIX PENCIL 243 
Then, since 6y = 0, we have at once 
b,=0, a@,=— by, a, = 5y,..., Gy —Op5, Ser, = 90, (10) 
so that [w™, —Aw™1, 2yw™-*,..., (—)-AKum-*, 0, 0,..., O]y = 0, (11) 
that is, [e*, —Aw*-1, uk... (—)FAK, 0, 0,..., Oly = 0, (12) 
and indeed 
[w*, —ApE-1, A2uk-2,..., (—)FAF, 0, 0,..., OJHA = 0, (13) 
since y was any column of HA. But this brings us to a minimal 
vector of order k < m, contrary to the hypothesis that m was mini- 
mal. Hence the assumption is false, and @ must be of the form 


6 = [u™, —Aw™—,..., (—)"r", 0, 0,..., OH. (14) 


3. Transformation of the pencil 
The matrix HA will now be studied. No element in its first row 
can contain pu, for the term in »”*! thus arising in the null vector 
6HA is impossible, having no other term in »”*! to cancel it. Hence 
the first row of HA is either a null row, in which case m = 0, 6 
being then [ 


or else this row contains at least one element aA. Let such an element 
be brought by interchange of columns to the leading position (1, 1), let 
the factor a be divided out from the first column, and let the A be 
used in elementary operations upon columns to clear any other 4’s 
from the first row. These operations are equivalent to post-multi- 
plication by a non-singular matrix K,, yielding at this stage HAK,, 
with first row [A, 0, 0,..., O]. 

Next consider the second row of HAK,. Under the A at (1, 1) there 
must be, for cancelling purposes, a p at (2,1) in an element which 
may be u+cA. There can be no other p» in row 2, for the Au” thus 
occurring in @HAK, could not be cancelled. If there is no A in this 
row either, then m = 1, for the minimal vector would then be 

[u+ca, —, 0, 0, 
with c = 0, because of the form of @ in (14). If m = 1, then HAK, 
is of the form, as far as its first two rows are concerned, 


k Tea fey | 
a ae ee ee 

Suppose, however, that there is a A in row 2. We bring it as before 
with unit coefficient to the position (2,2), and use it to clear other 
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\’s from the second row. We now have 

ae + i 

eas 

Under the A at (2,2) there must be, for cancelling purposes, a p at 
(3, 2) in an element of the form y-+-cA, and in the same way as before 
there can be no other p» in row 3. If there is no A either, we have 
A ; ae” 


Miu * 


and then m = 2, for 
[u(u-+ca), —A(u+cd), A?, 0, 0,..., 0JHAK, = 0, 
but again in such a case c = 0, in view of (14). 
Hence at this stage we have 


A 
HAK,=|# 2 - - - 


If there is a A in row 3, we bring it to (3,3) and carry out the same 
routine of elementary operations as before. 
The process must come to an end, and so, for a general value m, 


of m, we must have 
[x -, "4 


ee 


HAK — > ; = 7 


m, 


7 i as 
“Hy ni 
a a 








where L,, is a rectangular matrix of m,+1 rows and m, columns as 
shown above, the submatrix denoted by the asterisk being not 
necessarily null. 
If there is row-dependence in the submatrix A, we reduce it in 
similar fashion, obtaining 
Bay, ° 
- x «+h 


* * Pix 


My > M4, 
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and so we proceed until row-dependence in the A,’s is exhausted, 
At that stage A has been transformed into 








If there exists column-dependence in A, we reduce it in quite 
analogous, but transposed, fashion, obtaining+ 


, 
, 
ms 








r+s 
where A,.,, is either square and non-singular or else does not arise at 
all. In the former case it can be reduced to either of the canonical 
forms (the ‘classical’ and the ‘rational’) appropriate to non-singular 
pencils of matrices (3). 

As for the asterisked sub-matrices, such non-zero elements as they 
may contain can be reduced to zero by simple elementary operations, 
which have been described in (3). 

The matrix of type HAK that thus emerges is Kronecker’s canonical 
form for the singular matrix pencil (1), the invariants under equivalent 
non-singular transformations being the invariant factors (or if pre- 
ferred the elementary divisors) of the non-singular core A,,,, and the 
two sets of minimal indices m; and m; (2). 


+ L’ denotes the transposed matrix of L. 
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ON ANALYTIC FUNCTIONS REGULAR IN THE 
UNIT CIRCLE 
By MARY L. CARTWRIGHT (Cambridge) 


[Received 6 July 1933] 


1.1. Ler f(z) = f(re'®) u(r, 0)+iv(r, 0) 
be regular for |z| < 1, and let 
v(0) = 0. (1.11) 


Hardy and Littlewood* have proved that if, as r > 1, 
~ \ 1/A 
| u(r, @)\4 a : O((] r)- »), 


Myr, u) | : 


aT 





where a > 0, A > 0, then 

M)(r,v) = O((1—r)-*), 
asr—>1. If A > 1, the theorem reduces to M. Riesz’s theorem; it is 
then true for « = 0, and can be stated in terms of boundary func- 
tions. Making A > ©, we have the comparatively trivial result that, if 


u = O((1—r)-*), (1.12) 
where a > 0, as r> 1, then 
v = O((1—r)-*). (1.13) 


This result is not true for «a = 0. 
The problem before us here is to prove similar results depending 


on hypotheses on w,(r,@), where w, uif u > 0, and uw, = 0 other- 
wise. Carathéodory’s theorem shows that 
Dy 
fire!) —f(0) =" {4(R)—u(0)}, 
: R—r 


where r < R < 1, and A(R) = max u(R, A). It follows from (1.11) 


0<@<27 
that f(0) = u(0), and so, if w < 0 for r < 1, we have 


Dy 
M(r) = maak f(re®®)| < Rot 1 u(0)|. 
Choosing R = 3(1+,7), we have 
M(r) = O((1—r)*). (1.14) 


* G. H. Hardy and J. E. Littlewood, Journal fiir Math. 167 (1931), 405-23. 
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If u.(r,0) = O((1—r)-*), (1.15) 
this method gives M(r) = O((1—r)-*-*). (1.16) 
The function f(z) = —(1—z)-*, whose real part is negative, shows 
that (1.14) cannot be much improved; but a study of special func- 
tions leads one to suppose that (1.16) could be replaced by 


\ 


M(r) = O(l—r)") if a<1 (1.17) 
M(r) = O((l—r)-*) if a>1. (1.18) 
This is in fact true. 
It should be observed that, since 


| u(r, 0) d0 = u(0), (1.19) 


ry 


- 
M,(r,u,) and M,(r,u_), where u_ = u if wu < 0 and 0 otherwise, are 
of the same order. Hence theorems about M,(r, uw.) and M,(r,f) can 
be obtained immediately from the theorems of Hardy and Little- 
wood for conjugate functions. The theorems for M)(r,u,), where 
1 <A < , are more difficult. They will be considered in another 


» 


paper. 

1.2. The proof of (1.13) was shown to me by Professor Littlewood; 
and, as it has not yet been published, I shall begin by reproducing it. 
I next consider the case in which (1.15) is satisfied for « << 1. When 
this is so, we can afford to lose something in passing from (1.15) to 
(1.17), and so we can use a simple method. In proving (1.18) we have 
to use heavier methods involving the theorems of Phragmén and 
Lindelof. 

I should like to thank Professor Littlewood for suggesting the 
problem from which these investigations began and for much valuable 
advice. 

2. We shall prove the results in a more precise form than that 
stated in the introduction. In order to do this we have to fix the 
value of f at the origin; if the hypothesis is on |u|, it only remains 
to fix |v(0)| which has already been done by (1.11), but if, as in the 
later theorems, the hypothesis is on w,, we have to suppose that 
u(Q) > 0, or some such hypothesis. In what follows A is an absolute 
constant; B, C, D, K may depend on a parameter ¢, say, in which 
case we write B(f). A, B, C, D preserve their identity throughout 
a theorem and its proof, but K need not be the same in each place. 
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THEOREM 1. Jf |u(r, @)| < A(1—r)-* 
where « > 0, for r < 1, then 
lw(r, @)| < K(a)A(1—r)~° 
forr< 1, 
Poisson’s formula is 
u(R, d)(R?2—r?) dd 


u(r, 4) = Qa aE R?+-r?2—2Rr cos(¢—8)’ 


where r < R < 1: and so we have 


dv leu ais R, b)( R?—r?)sin(d—6@) dd 
{ R?-+-r?— 2 Rr cos(d—@)}” ; 


or 6 Tv 
7 


We may without loss of generality suppose that 6 = 0, and we have 
l [ sin | \dd 
27 J {(R—r)?+4Rrsin?3¢}* 


« 


A(1— R)-*(R2—r?)- 


Putting R = 3(1+,7), we have 


| a9] db ¢ 
"| < KA(1—r) “fe oe 


ori 


pdp — 


Integrating, we have 
lu(7, 0)| < K(a)A(1—r)-°, 
which is the required result. 

3.1. In the proofs of the two following theorems we subtract a 
function whose real part is positive so as to obtain a function whose 
real part is negative in a certain region. We then apply Carathéo- 
dory’s theorem in that region if possible; if not, we map the region 
conformally on a more convenient one, and then apply Carathéo- 
dory’s theorem. 

THEOREM 2. Jf u(0) > 0, and if 

u(r, 0) < A(1—r)-4, 
where «a <1, forr < 1, then 
M(r) < K(a)A(1—r)" 
forr <1. 
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We first observe that, applying Carathéodory’s theorem to f(z) in 
the circle |z} < R < 1, we have 


2r 
Iflre®)—f(0)| <p" {A(R)—u(0)}. 
Putting R = }(1+r), and remembering that 1(0) = 0, u(0) > 0, we 
have | f(re')| < 4A(1—r)-*2+.|f(0)| < 4A(1—r)-"1-4.4. (3.11) 
It is sufficient to consider a typical point, z = re%, and show that 
f(re’®)| < K(«)(1—r)“, where K(a) is independent of 4). There is 
no loss of generality in supposing 0, = 0. Let 
fil) = uy +tv, = f(e)—KA(1—z), 
where a < 8 < 1, and consider the value of u, inside the region 
\1—z| < 24(1—r)™B, 0 < |1—z| < Ry. (3.12) 
Writing 1—z = pe’, we have 
R{(1—z)-*} = p-Pcos Bb > p-Feos $8 
inside the unit circle, and so 
R{(1—z)-?} > 2-48(1—r)-*cos $87 
in the region defined by (3.12). Choosing K > 2!sec 87, we have 
u, < u—KAR(1—z)-8 
< A(1—r)-*—A(1—r)-* = 0 (3.13) 
in (3.12). 
If « < 4, we can choose B > 2a, and then it is easy to see that 
the region (3.12) includes the circle |z—4| < 4, provided that Ry > 1. 
For inside this circle r < cos@, and so 


|l—z|? = 1—2rcos0+r? < 1—?r? < 2(1—7). 


In this case, applying Carathéodory’s theorem to f,(z) in the circle 
z—4| < 4, we have 

u,()| 

fry < Os pani 
and it follows from . that 


\1(4)| < |f)| < (2-*+84+1)A+K()A 28 < K(a,B)A. 


Hence Lflr)| < I f,(r) K(B)A|1—re®®|- = 
< nee 1 = K(a)A(1—r)** 


which is the required result. 
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If x > $, we put 
1—z = (1—w){1—(1—-w)y", (3.14) 
where 0 < y < 1, 1—z = pe’?, 1—w = ae’#. Then for appropriate 
Ry. R), y, the region defined by 


b|<4x, o< RP 


1 Ri (3.15) 
is included in that defined by (3.12), viz. 
p < 2(1—r)¥B, p< Rp. 


oYsin yi - 
sO 


For, as p> 0, 0 ~ p, and ¢ ys— yx, where tan x 


1—oYcos 


x ~ p’sin yy. 
Hence || < 47 is practically 
—Jo < $-+prsin yp < de. 
-47, so is &. It follows that 


9 


If d is near }z, so is #, and if ¢ is near 
the region (3.15) is included in 
dar— |b| > SpYsin hyn, px Ry (3.16) 
for sufficiently small Rj. But, since @ < 29, as z > 1, we have 
kar—|b| < 2cosd = 2(1—rcos @)|1—z|- < 2(1—r+4r6?)|1—z|"1 
2(1—r)p-1-+ 2rp (3.17) 
for sufficiently small p. It follows that the region (3.16) is included in 
Sp¥sin ym < 2(1—r)p-!+ 2rp, p< R,, 
i.e. p’*+1— 2p?cosec Sym < 4(1—r)cosec $yz, p< R,. 
Since y < 1, this region is included in 
py aoe 
and, if y+1 < B/a«, this region is-included in (3.12) for sufficiently 
small R). 
Now we apply Carathéodory’s theorem to 
O(w) = f,{w+(—w)r*}} 
in the circle jw—1+4Rj| < }R;. It follows from (3.13) that 
R{®(w)} < 0 in this circle, and so 
2R{P(1+ LR) O(1-+4R))). 


2 
] Ww sii. 


8(1—r)cosec yz, p< Ry, 


D(|w!) 


Since R; depends only on a, f, y, it follows from (3.11) that 
@(1+-3Rj)| < K(a,B,y)A = K(a, B)A; 
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and, since |l—w| ~ |1—z|, we have 
filr)| < K(a,B)A(I—r)4 

from which the required result follows as before. 

3.2. For a > 1, the method gives the following result, which will 
be superseded by Theorem III: Jf u(0) > 0, and of 

u(r,@) < A(1—r)-*% 
where « > 1, then, given any « > 0, we have 
M(r) < A(1—r)-*-€ 

provided that re) <r <1. 

We have 

R{(1—z)-*} > cos aB|1—z|-* > cos aB(cos B)*(1—r)-* 
for |jarg(1—z)| < B < 4m/a, |1—z| < R, = R,(f). Hence 
aa )} = R{f(z)—A sec af(sec B)*(1—z)-*} < 0 
in that region. Putting 1—z = (1—w)?f!7, we have 
R{O(w)} = R{f,(1—(1—w)*#l*)} < 0 
in the cirele }w—1+4R/| < }R; for sufficiently small Rj. Applying 
Carathéodory’s theorem to ®(w), we have 
2R{M(1+}Ry)} 
i—| |w | 
But Rj depends only on f, and (3.11) still holds. Hence 
(\w))| < K(B)A(L— |), 


M(\w))| < 4+@(1+4R’). 


and s¢ 
wii < |f,(r)|+A sec aB(sec B)*(1—r)- 
< K(p)A(1—r)-Al= 
for R’ <r <1. Given any « > 0, we can choose 8 so near to $2/a« 
the ut Bn <a+4e, and then Rj so near 1 that (l1—r)-* > K(f) for 
< 1, and we have the required result. 
1. It remains to prove the more precise results for a > 1. 


TueoreM III. Jf u(0) > 0, and if 


u(r, 0) < A(1—r)-* 
where « > 1, then 


u(r, 0) > —K(a«)A(1—r)-* 
and so M(r) < K(a)A(1—r)-*. 


This theorem is an immediate corollary of the following lemma. 
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Lemma 1. Suppose that F(z) is regular and 
log| F(re)| < Br* 
for r >1, \0| <B, where B > 3n/a. Suppose also that F(z) has no 
zeros in that region.* Then, given any 5 > 0, we can choose K(8) so that 
log|F(re®®)| > —K(8)Br* 
for \6| < B—8, r > 1. 

For if we consider the typical point z = 1, and put 1—z = 1/w, 

we find that, since a > 1, 

F(w) = exp{f(1—1/w)} 
satisfies the hypotheses of the lemma for some f such that }2/a < f 
< 4m with / = secf8, and B = K(f)A. Hence, putting 5 = $8, we 
have 

u = log|F(w)| > —K(B)A|\w|* 
= —K(a)A|1—z|-* > —K(a)A(1—r)-* 

on the real axis, which is the required result. 

4.2. The lemma depends chiefly on the theorems of Phragmén and 
Lindeléf.+ If B + 0 as 1 — o, the lemma reduces to a simplified form 
of a theorem which I have given elsewhere, and the proof for B > 0 
follows the same lines. The fundamental theorem of Phragmén and 


Lindeléf on which nearly all the others depend may be stated as 


follows: 
THEOREM OF PHRAGMEN AND LINDELOF. Suppose that F(z) is 
regular anc 
gular and log|F(re!)| < Br 
for |0| < B < 4n/a, r >I, and let 
— log | F(re*®)| 
A(6) = lim —2 sail C 


na 
Suppose further that a 
(0) < Ccos0a+Dsin bx 
for 6 = +8; then, given any « > 0, we have 
log| F(re’®)| < {H(0)+e}r* < (C'cos 0a+D sin 0a+«)r* 
for \6| < B and every r > r,(e). 


* If F(z) is not identically zero, but has some zeros in the sector, it is still: 
possible to obtain a result of this kind which is valid except in certain circles 
near the zeros. Compare M. L. Cartwright Proc. London Math. Soc. (unpub- 
lished). 

+ E. Phragmén and E. Lindeléf, Acta Math. 31 (1908), 386. 

t M. L. Cartwright, Proc. London Math. Soc. (unpublished). 
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Proof of Lemma 1. We may, without loss of generality, suppose 
that 8 < 7/«. For if not, we can obtain the required result by dividing 
the angle |@| < 8 into a number of angles each less than 27/a and 
applying the result for 8 < 7/« in each of them. We have H(#) < B 
Suppose first that H(y) = —oo, where |y| < 8. Then, applying the 
above theorem in the angles y, < 6 < y, y < 6 < ye, where —B < yy, 
Y—-V < 7/0, ye—y < m/a, ye <B, we see that H(?) = —oo for 
y¥, < 9< yp. And, repeating the argument if necessary so as to cover 
the whole angle, we have H(@) = —oo for |@| < B. Since B > z/a, it 
follows* that F(z) = 0, which is contrary to our hypotheses. Hence 
H(@) is finite for |@| < B. 
In this case, since }7/a < B < pe we havet 
H(0) > M{H(B)+H(—A)}sec Ba. 
For if H(0) = }{H(8) His B)}sec Ba—, (4.21) 
where 0 < { < ©, we hav ve 
H(@)< <iteg- ae B)}sec Ba cos Aa+- 
+ }{H(8)—H(—£)}cosec Basin @a+{ cosec Basin(@—B)a (4.22) 
for 6 = 0 and 6 = Hence (4.22) holds for 0 < 6 < B. Similariy, 
H(6) < }{H(8)+A(—B8)}sec Ba cos 0a+ 


+ }{H(B)—H(—B)}cosec Bau sin @a—{ cosec Ba sin(B+-0)a (4.23) 
for —B <6 <0. But it follows from (4.22) and (4.23) that 
1 


H(0) < 4{H(8)+H(—B)}sec Ba cos 0a+- 
+ }{H(8)—H(—B8)}cosec Basin 0a— 
—{cosec Busin(B—y)asec yacosOa, (4.24) 
where 0 < y < 4a/a, for 0= +y. Hence (4.24) holds for |@| < y 
which contradicts (4.21) if y is chosen equal to B—}7/a. 
Consider the function 
O(z) = F(z)exp{(C+iD+Q)z%exp(—érs), 
where ( > 0, € > 0, 
C = —}{H(B,)+H(—B,)}sec Ba, 
D = }{H(B,)—H(—B,)}eosec B, a, ; 
and 42/a < B, < B. Now, given any e, > 0, we can choose ry = 1(€;) 
so that log | F(r re’) )| < {H(6) )t+e,}ro 


* See Pélya, and Szegé, Aufgaben und Lehrsdtze, i (1925), 148, No. 327. 
+ Compare G. Pélya, Math. Zeitschrift, 29 (1929), 574. 
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for |@| < B,, r > ro(e,); and so, if «, < a. x, we have 


log |®(re'*:)| < [H( ‘o wr {H(B —B,)} 

—4{H(8,)—H(—8,)}-4 pata ipo < (e,+fcos B,«)r* < 0 
for r > ry, and pelle. log |®(re-‘*:)| < 0 for r > 79. We choose 
so that 
for |0| < B,. 

It follows from the definition of H(@) that, given any «, > 0, we 
can find a sequence R,, R, R,, > «© such that 

log|F(R,,,)| > {H(0)—e)}R*, > [4{H(8,)+H(—B,)}sec B,a—e.|R 
Hence log |®(R,,,)| > (C—e,) R* —&rf > 0, 
provided that «, << { and m > m,(e;, €,€). We apply the maximum 
modulus principle to ®(z) in the region ry <r < R, R> R,,, 

Since |®(z)| < 1 for |z Y), |0| < B,, and for 6 = +,, 

R, while |®(R,,)| > 1, there must be a point, Re’®, on the 

R, \0| < B,, at which 

O( Re'?)| > 1. 
Taking this point as centre and applying ( ‘arathéodory’s theorem* 
to log ®(z) in the circle |z— Re'®| < Rsin(S—f,), we have 
log |®(z)| > —K(B—B,)BR* 
for |z— Re’®| < }Rsin(8—f,). That is, 
log |®( Re'*)| > —K(B—B,)BR* 
for |@—O| < 4(B—8,). 
Let ©, and ©, be the two extreme points of the arc 
9—O| < 3(8B—B,) 


which lie inside the are |6| < B,. One of them may coincide with 


log |®(r, e!”)| < 0 


B, or —f,, or both may lie inside |6| < f,. Using (4.25), we apply 
Carathéodory’s theorem to ®(z) in the circles 
z— Re'®:| < Rsin(B—B,), z— Re'®:| < Rsin(B—A,), 
and it follows that (4.25) holds for 
6—O,| < 4(B—B,), 6—O,| < 3(B—f,). 

Hence (4.25) holds for that part of the are |@—O| < B—f, which 
lies inside |@| < B,. Repeating the process not more than 48,/(8—,) 

* Perhaps it would be more correct to say that we apply Valiron’s theorem 


to M(z), see G. Valiron, Lectures on the General Theory of Integral Functions, 
Cambridge (1923), 89. 
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times, we find that (4.25) holds for || < £,, R > R,,,, which is 
equivalent to the desired result. 

4.3. For « = 1, we have the following theorem: 
THEOREM IV. Jf u(0) > 0, and if 
u(r,0) < A(1—r)" 


then Mir) < KA(—r)-(log : : (4.31) 
l—r 


The number 2 in (4.31) cannot be replaced by any smaller number, 
as may be seen by considering the function 


fz) =—- Z (log +.) 


1—z 


Putting 1—z = pe?, we have 
eeu — 8 (log ) +¢2\4 apa, ! 
p |  s 
If cosd > xf{log(1/p)}, then wu < 0, and if 
0 < cos¢d < aflog(1/p)}, 
then 1—r < azpf{log(1/p)}-. 
Hence u < (1—r)" 
for r < 1. But obviously 
M(r) ~ 2 (los . J 
l—r 
as f > 0. 
The theorem is a corollary of the following lemma, which corre- 
sponds to Lemma I. 


Lemma 2. Suppose that F(z) is regular and 
log| F(re"®)| < Br 
for r > 1, \0| < 4a. Suppose also that F(z) has no zeros in that region. 
Then, given any 5 > 0, we can choose K(8) so that 
log|F(re®)| > —K(8)Brlogr 
for \0| < 4n—6, r > Rp. 
For considering the typical point z = 1, we have 
2 1 ] 
u< KA i 12 '°8 ia] 
in the region 


"A g(a eons 
l—-r>K\l |\lo8: a’ 
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Putting 1—z = (1—w)log1/(1—w), 1—z = pe'?, 1—w = ae'¥, we 
find that the region (4.32) includes the region 

lb| < 4z, o< R, (4.33) 
provided that K and Rj are suitably chosen. For p ~ oflog(1/c)} as 


p—> 1, and 1 | 1\-1) 
= == J 1 rc~ Té ] PF ne , 
o—y arg | og — tan ( og ; | 


= 4 
~ —v(lox) ~ —v(loe:) . 
0 p} 
Hence (4.33) is included in 
PS Ro, 
1.e. 1 \p| —1 e<R,; 
and so it follows from (3.17) that (4.33) is included in 
47(log 1/p)-! < 2(1—r)p-1+-2rp, p< R, 
which is included in (4.32) for K > 32, Ry = R,(K). 
Hence in the region (4.33) 
bt Che 
J |l—w 


' 


R{@(w)} = R{f(1-C—wlog, : 


— 


] 


It follows that F(w’) = exp((1— =) 
Ww 


satisfies the hypotheses of the lemma for some B, and we have 
u = R{@(w)} = log| F(w’)| > —KAw’ log w’ 

= ~ Kien — 6) ee 

1 *|l—w |] 


|l—w 


for jarg(1—z)| < 47—8, |l—z| < Rj. Hence 
. ] : } 
u(r,0) > —KA log : 
l—r l—r 
and it follows that 
] 1 \2 
(r,0)| << KA lc 
0(7,0)| < KA>* (logs) 
by the method used in Theorem I. 


4.4. Lemma 2 depends on a simple special case of Carleman’s 
formula* which may be stated as follows: Suppose that F(z) 1s regular 
and has no zeros for \z| >1, |argz| < $a. Then 


* See G. Pélya arid G. Szegié, Aufgaben und Lehrsdtze, i (1925), 120, No. 178. 








ON ANALYTIC FUNCTIONS REGULAR IN THE UNIT CIRCLE 257 
4a 
ee a 
aR | log | F(Re'®) cos @ d@ + 
har 


R 
lf toe! Ftér\| loo F(— arnt — 2 
+ oF | {log | F (ir) |+-log| F(—ir) i(a— 7a) af 
i 


bounded as R > co. 
Proof of Lemma 2. Since log|F(+ir)| < Br for r > l, we have 


1S 


aH 


R R 
DD : a beeen ae : cs B fd B 
= | flog | F(ir)+log F(—in) (5-75) dr < = | af < 7 log 5. 
i 
Hence (1/7R) [ log| F( Re) \cos 0 d0 > —(B/z)log R+K,. (4.41) 
—_ir 
It follows that, given any « > 0, we can choose K(e) so that 
log| F( Re'®)| > —K(e)BRlog R, (4.42) 


except in a set of 0 of measure less than ¢, for |@| < $x, R >I. For 
if log| F( Re®®)| < —K BRlog R in a set of 6 of measure ec, then 


(1/7R) | log|F( Re’) |cos 0 d0 < —(KB/z)log R | sin 0 dé 
E 0 
< —(KB/z)cos clog R, (4.43) 
where E denotes the set of 6 of measure «. The integral taken over 
the complementary set is certainly less than 
(B/z) { cos 0 d0 = 2B/n, (4.44) 
= 
and so, if Keose > 1, (4.43) and (4.44) give a contradiction of 
(4.41) for R > R,(e, K). Now we apply Carathéodory’s theorem to 
@(z) = {log F(z)}/Rlog R, 
in the circle |z— Re’®:| < Rsin(}z—O,), where ©, belongs to the set 
for which (4.42) holds and |0,| < «. Given any k > 1, we have 
log |®(z)| > —K(e,k)B 
for |z—Re'*| < Rk-sin(}7—@,); and, given any 6 > 0, we can 
choose « and k so that this circle includes the are |#| < $7—8, 
Rsin6d for all R > Ry. Hence 
log| F(Re’®)| > —K(e, k)BRsec8 log( R sec) 
> —K(d)BRlog R 


for R > Ry, which is the required result. 
3695.4 Ss 
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THE ENERGETICS OF NON-STEADY STATES, WITH 
APPLICATION TO CEPHEID VARIATION 


By E. A. MILNE (Ozford) 
[Received 13 October 1933] 


1. THE following paper develops a new method for discussing non- 
steady states of thermal systems. A thermal system is in a ‘steady 
state’ when the temperature at every point, not necessarily uniform, 
remains constant in time. In that case the loss of energy in each 
small element in any interval of time due to any causes must be 
balanced by an equal gain due to other causes. Integrating over the 
whole system, we have the result that the rate of loss of energy, at 
the surface, or at internal bounding surfaces delimiting ‘foreign’ 
bodies not parts of the given system, must equal the rate of gain of 
energy from given sources. The simplest example is an electric fire, 
the source of energy here being the applied voltage; another example 
is a star in a steady state, where the radiation to space is balanced 
by internal sources of gravitational or sub-atomic energy. The pre- 
sent paper sets out a general method of solving the question of what 
happens when a system in a steady state is disturbed, it being 
assumed that the supply of energy is continued at the same constant 
rate as before. The theory should apply to Cepheid variable. stars. 
The results show that according to circumstances the system may 
or may not return to the steady state; it may return to it by ‘damped 
oscillations’, or in a ‘fully damped’ manner, or it may diverge from 
it, in a systematic direction or by oscillations of increasing amplitude. 
But the most interesting result is that the disturbed system may 
under certain conditions oscillate with a definite period, and that in 
this case the oscillations are definitely of non-harmonic type, pos- 
sessing maxima of energy-output which are more pronounced than 
the minima. This asymmetry is a well-known characteristic of the 
light-curves of Cepheids, which thus receive a measure of explanation. 

The theory is concerned solely with energy considerations. Just 
as in dynamical problems of one degree of freedom, in which the 
motion can be obtained from the constancy of total energy without 
consideration of reactions, external or internal, so the present theory 
fixes attention on the constancy of rate of input of energy and 


v 
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determines the nature of the output. All consideration of internal 
circumstances such as pressures, temperatures, etc., is avoided. The 
systems examined differ from conservative mechanical systems in 
that the total energy stored in the system is not constant in time 
but fluctuates. 

The analysis hinges on a functional equation of the type 


K$(t) + b{t+b+4()} = 0, 

and the paper has a mathematical interest of its own, apart from its 
physical applications, in the classification of solutions 4(¢) of this 
equation according to the values of the constants « and b. The con- 
stant « plays the part of a damping parameter; its value determines 
whether the system is oscillatory or non-oscillatory, stable or un- 
stable, and if oscillatory then whether it is periodic. The constant 
b determines the period in the case of strictly periodic disturbances, 
and the quasi-period in the case of general oscillatory disturbances. 
The use of the functional equation may be considered as alternative 
to that of exponential functions and Fourier analysis, which play 
here no part. It is remarkable that the above simple functional equa- 
tion itself generates the classification of disturbances classically 
treated by the Lagrangian method of using exponentials; its advan- 
tage is that it generates periodic functions directly, and not as the 
sum of a number of harmonics. 


2. Consider a system which generates energy, or to which energy 
is supplied, at the constant rate H. Choose @ zero of time. Then at 
time ¢ after this zero, the total amount of energy generated or sup- 
plied is Ht. 


3. Let f(t) be the total amount of energy lost by the system in the 
same interval (0,#). It may have been emitted to space, as in the case 
of a star; or converted into mechanical energy and then back into 
heat by friction with bounding surfaces, as in the case of ‘Griffiths’s 
Heat Engine’ ,* or in other ways. If the state of the system is steady, 

f'(t) = #, and the total energy imprisoned in the system remains 
constant. By integrating we have 


S(t) = #t. (1) 
* E. H. Griffiths, Thermal Measurement of Energy (Cambridge, 1901), 50. 


I am indebted to Professor H. F. Newall for this reference, and for introducing 
me to this example of the general theory. 
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In a non-steady state f(t) will not be in general equal to Ht. The 
energy imprisoned in the system, say Q(¢), will not be constant in 
time, and by counting up gains and losses we have 


Q(t) = Ft—f(t)+Q,, (2) 
where 2, is the energy imprisoned at ¢ = 0. 


4. As time goes on, in general f(t) steadily increases. In the case 
of a star, this is obvious, as the differential coefficient f’(t) is just 
L(t), the absolute bolometric luminosity, which is always positive. 
The following arguments do not depend on the assumption that f(é) 
steadily increases, but it is clear that in the neighbourhood of any 
stable steady state f(t) roughly keeps pace with Ht. Corresponding 
to any epoch ¢, there will be an epoch ¢+7', where 7’ is a function 
of ¢ such that f(¢+ 7), the energy emitted up to time t+ 7’, differs 
from Ht only by an amount depending on the initial conditions. We 


write therefore ft+T) B(t-+e), (3) 


where c is a constant, depending on the conditions at t= 0. It is 
clear in fact that the early behaviour of f(t), near t = 0, may be quite 
independent of the rate of supply; for the system may be started in 
a state far removed from the steady state corresponding to #, and 
its initial radiation to space will depend on this initial state. The 
physical meaning of ¢ is seen by putting ¢ = 0 in (3); we get 


Sf T)) Ec, 


where 7), is the value of 7 at t= 0. Thus Ec is the emission of 
energy between? = Oandt = 7). Roughly speaking 7) is a measure 
of the time that must elapse before the supply influences the emission, 
but this is only a crude interpretation. The definition of 7’ is simply 
(3). 7’ may be called the time-lag.* The idea of a time-lag in Cepheid 
variation between the generation of energy in the interior and its 
emergence at the surface is not new; it has been specifically men- 
tioned by Eddingtont and by Russell, Dugan, and Stewart.{ Edding- 
ton definitely contemplates, moreover, the possibility of H remaining 
constant during the variation. 


* *Time-lag’ occurs in the theory of elastic fatigue, and has been treated 
by Volterra by the method of ‘functionals’. 
+ Internal Constitution of the Stars (Cambridge, 1926), 202. 


t Astronomy (Grim, 1927), 2, 767. 











THE ENERGETICS OF NON-STEADY STATES 261 
5. As a particular case, when the state is steady we have, by (1), 
ft+T) = Et+7), 
so that comparison with (3) gives 
T=€c 

and 7 is constant, as is physically obvious. For this case we write 
T =F 

6. The reader must recognize that the definition of 7 is simply 
equation (3), and not seek to impose any more intuitive meaning on 
T. If we could regard elements of energy as maintaining an indepen- 
dent identity during their career from generation or supply in the 
interior to emission at the surface, 7’ could be regarded as the mean 
time of transit. But this is not legitimate. We could, of course, 
attempt to calculate the delay between a sudden temporary increase 
in 2 and the peak of the consequent increase in surface emission in 
an otherwise steady state, but this is not exactly 7' as defined. As 
it is, I propose to regard 7' as having a physical meaning, and as 
being a function of the state of the system at time ¢. Equation (3), 
then, definitely assures that sooner or later all the energy supplied 
to the system by time ¢ is subsequently disposed of. I proceed now 
with the consequences of definition (3) and the assumption that 7’ 
depends on the state at time ¢. 


7. In the general case of a non-steady state, the time-lag 7’ and 
the imprisoned energy Q are both functions of ¢, and we may there- 
fore represent 7' as a function of Q. Jt is assumed that T is a one- 
valued function of Q. In the case of periodic oscillations, for example, 


in which Q passes from Q,,;,, to Qyax and back to Q,,;, in a period, 
the assumption means that 7' is supposed to take the same value 
at the two epochs at which Q takes any given value lying between 
QO. im ANd Q,,a%- It does not seem possible to offer any a priori physical 
justification for this assumption, though the assumption itself appears 
not unreasonable simply as a hypothesis whose consequences are to 
be investigated. Actually the consequences are of sufficient interest 
to afford some a posteriori justification for making the assumption. 
Further, the single-valuedness of 7' as a function of Q is readily veri- 


fied in the later analysis describing the resulting physical situation. 


max min 


* For further discussion of the meaning of 7,, see § 12 below. 
+ I am indebted to Professors D. R. Hartell and R. H. Fowler for their 
discussions on this point. 
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We do not, however, assume that the system is stable in its steady 
state, but consider the most general small departures from the 
steady state and investigate their subsequent behaviour. Write, then, 

T = F(Q) 
= F{Q,+(Q—Q,)}, (4) 
where Q, is the energy imprisoned in the steady. state. Since Q—Q, 
is a small quantity, we have approximately 
T = F(Q,)+(Q—Q,)F’(Q,). 
Substituting for Q its expression (2), we have 
T = F(Q,)+{#t—f()+Q,—Q,} F’(Q)). 

Write this as 


i Kee . ‘ 
T ‘ ie git )—E;, (5) 
where y and « are constants. Their values are given by 
y = FQ,)+(Q,—Q,) F’(Q,), 
or approximately , 
PI =" > I (Q25) 
1) (6) 
and by x = —EF'(Q,) 
(aT * 
, (aa), re “) 
8. We now substitute from (5) in (3), obtaining 
a . . 
lt+y+GU0— BH |= Be+e) (8) 
Define a function ¢(¢) by the relation 
: . t 
fi) = E (« 14% ), (9) 
K 
a being a constant to be chosen later. Then (8) reduces to 
l 
a(l+x«)+y+¢(t)+-¢[t+y+Ka +-(t)} Cc. (10) 
K 
Choose a so that tn. (11) 
l+-K 
Then y+Ka Mos (12) 
1+-K 
We put a ah, (12’) 
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and (10) then becomes 





| « $O+${l+b+o)} = 0 |. (13) 
This is the fundamental equation which ¢ must satisfy. In terms of 
¢, (5) gives for T' T — b+¢(t). (14) 


Thus the physical meanings of ¢ are: (a) that it determines the 
fluctuation of 7’ with time; (b) that it determines the fluctuation 
with time of the difference between f(t) and its secular term Et. If 
we write f’(t) = L(t) (adopting the usual astrophysical notation), we 


have from (9) ’ 
L(t) = w+ ©). (15) 
K 
Writing bL(t) = Lit\)h—E 
and E= iL,, 
where L, is the steady-state luminosity, (15) gives 
ski) _ #() aay 
L, K 


4(t) has the dimensions of a time, and « is dimensionless. 

9. From physical considerations ¢(¢) is a continuous differentiable 
function of ¢ in any finite range of ¢ for ¢ > 0. The behaviour of ¢ 
depends on x, and we classify the cases as follows: 


Case 1. « = 1. Equation (13) becomes 


$(t)+bt+b+4()} = 0. (17) 
teplacing ¢ by t+b+¢(t), we find 
$(t+-2b) = g(t). (18) 


Hence ¢ is periodic, of period P = 2b. Putting « = | in (12’), the 
value of P is y+c, or T,+7;, so that b is essentially positive.* We 
shall give yet another physical meaning to 6 later.+ 

We next note that ¢(¢) cannot take the value —b. For if, say, 
A(t’) = —b, then by (17) —2b = 0, which is a contradiction.{ Since 


* It is readily verified that the main result about to be obtained, namely, 
that the interval between minimum and next succeeding associated maximum 
of ¢ is less than the interval between the maximum and the close of the period 
(next associated minimum) holds good for functions ¢ satisfying (17) even if 
b is a negative number. 

+ We show that b is the time-average of 7’. 

* We exclude the special case b = 0 for the time being. 


+ 
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¢ is continuous and periodic, there must be a range (—b—y, —b+-8) 
(yn, 5 > 0) inside which ¢ cannot lie; again, since it is continuous, it 
must lie wholly on one side of this range or wholly on the other. 


Hence either 4(t) < —b—y or d(t) > —b+68. If g(t) < —b—7, let 
it take its upper bound —b—7y for t=t’. Then by (17), putting 
t=, (t'—n)=—b+y, contradicting ¢(t) < —b—y. Hence 
d(t) > —b+6. The upper bound of ¢ cannot exceed b—8. For if, 
say, (t’) M > b—8, putting ¢t = ?’ in (17) we have 

d(t’+b+M) M < —b-+4, 
contradicting 6 > —b+6. The upper bound of ¢ is, in fact, precisely 


b—8. For if the lower bound —b-+-3 is attained for t = t,, then, by 
(17), d(t, +8) b—6, so that the value b—6 is attained for ¢ = #t,+6. 


It follows that —b+8 < dit) < b—8, 


whence 0<s8<b. 
© 


Should 6 = b, the solution reduces to ¢(t) = 0, and the state is 
steady. 

Thus, provided § 4 b, every minimum —b-+ of ¢ is followed by 
an associated maximum at time 6 later, and since b = $P, 6 is less 
than half the period. Accordingly the interval from minimum to next 
succeeding associated maximum is less than half the period. This asso- 
ciated maximum is easily shown to be followed by an associated 
minimum at the end of the period (reckoned from the first selected 
minimum), and accordingly the interval from maximum to next 
succeeding associated minimum is greater than half the period. There 
may be any number of minima —/+6 and maxima b—6 in the 
period, but they may be associated in pairs in the manner 
indicated. 

The graphs of a typical function ¢ and its derivative are shown 
below (Fig. 1); ¢’ is here the proportional luminosity fluctuation, 
dL(t)/L,, whilst —d¢, apart from a constant, is a measure of the 
imprisoned energy, since 

Q(t) E{a-+-¢(t)}. 
The steady-state luminosity L, is also the mean value of L(t). For 
LP +P 
; 1+ 
Lit) dt = 5 


ty t, 


> E{1+-¢4'(t)} dt = Ly. 
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We have the further quantitative relations 
ty 


fu 


7 t= $+8)—$(h) = 26-8) = P—28, 
a | 





ae Y : 
p(t) cr> 7 “ wid 
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Hence in the portion of the period, from zero to zero, in which 6L(t) 
is on the average positive, 
— i 
dL(t) _ P 3%. (19) 
Ll, ) 
and in the portion of the period, from zero to zero, in which 6 L(t) is 
on the whole negative, 





—T 
+ 
mn 





Fie. 1 


rs (20) 


\S.L(t) 
L,  P-3 
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These are capable of observational tests on Cepheid variables. In 
terms of apparent bolometric stellar magnitudes m we have 


L, = Alo-4m, LE = A10-04m, 


where A is a constant, and m, is defined by 
L+P 
10-0°4 ae ; 10- 0-4m dt. 


th 
Relations (19) and (20) then give 
4,48 


‘ 1 4 
a [ 1Q-0-4(m 5 (19’) 


P 
190-0 4m m,) dt = 3 = (20’) 
P—5 


It follows at once from these relations that the mean amplitude 
of the are of the curve 6L(t) which rests on the shorter base A’ B’ 
(see Fig. 1) is greater than the mean amplitude of the are which rests 
on the longer base B’C’. This refers to curves on an absolute energy- 
scale. It is well known that the light-curves of Cepheids on a magni- 
tude scale usually show more pronounced maxima than minima, in 
general agreement with the present theoretical predictions. It may 
be mentioned that this feature is to be distinguished from another 
feature of Cepheid light-curves, namely, the steeper rise to maximum 
than descent from maximum, which is not accounted for by the 
present theory. 

Let t, be a zero of d(t). Then by (17), putting ¢t = i), we have 
d(t)+b) = 0. Thus the interval between two associated zeros of ¢ is 
half the period. It follows that 

to 4,48 
[ L(t) cf [ 5 L(t) 


ad | ed L, 


dt = 4P—1 (21) 


t 
| ee | (22) 


ty 4P 


These are also capable of observational test; tf) is the abscissa of the 
vertical which divides the lobe A’B’ of the energy-curve into two 
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equal areas; this abscissa should differ by }P from that of the vertical 
which divides the lobe B’C’ into equal areas. 

The mean value of ¢(¢) over a period is zero. For 


t,4.2b +8 t£2b 
| ¢oa= f[ gmat | goae 
ty ty 4,48 
4,46 t,42b 
=— f[ gtt+o+go}a+ | god. 
i, t,18 
In the first of the integrals on the right-hand side, change the variable 
from ¢ to 6, where 


t+b+d(t)=0,  {1+¢'(t)} dt = dé. 


By differentiating the functional equation for ¢, we find at once 


(1+4'(}{1+4'(6)} = 1. 


This relation shows that 1+¢’(t) can only change sign if 1+¢'(@) 
becomes infinite, which is excluded on physical grounds. Hence 
1+-¢4'(t) is always positive or always negative, and since it is unity 
at a maximum or minimum of ¢, it is always positive. The trans- 
formation of variables is therefore legitimate.* The limits for @ are 
?, +6 and ¢,+2b. Hence 

t,4-2b t£2b t-£2b 

| ga =— | gO1+¢(O}d0+ | ge)de 

t, +8 t,+8 


2 O=t,+2b 
6=t,+8 


since ¢ takes numerically equal values at the limits. It then follows 
from (14) that 7’ = b, ie. the period P is twice the mean value of 
the time-lag. There is a sort of resonance between the oscillations 
and the time-lag. 

The functional equation in effect maps the are BC of Fig. 1 on the 


are AB. For astvariesfrom t, to 4+8, 


d(t) varies from —b-+6 to b—6, 
t+b+¢(t) varies from ¢,+5 to t,+26, 
d{t+b+¢4(t)} varies from b—8S to —b+6. 


* T owe this point to Professor Titchmarsh. 
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In t, < t < t,+8, we may choose for ¢(¢) any function satisfying the 
conditions $i} = Sas, $(t, +8) = b—8, 
$'(t,) = $'(4,+8) = 9, 
and otherwise arbitrary; in t,+65 < t < t,+2b, ¢(t) is then defined by 


ot) = —9(") (A> 


where ¢’ is the appropriate root of 
t’+b- b(t’) t. 


Thus discussion of the energetics of the system, though fixing certain 


periodicity properties, leaves a considerable degree of adjustment 
open for the details of the particular system considered. 

It should be noted that by applying the inequality —b < ¢4(t) < b 
to the particular case t = 0, at which, by (9), 6(0) = —a, we have 

(y+) - (c—y) < ye, 

which requires y > 0, ¢ > 0. These conditions are consistent with 
the physical meanings of c and y, namely, c = 7, y = T. They 
also show why it was necessary to introduce c. 

The case « = | is the most interesting one physically. We proceed 
now to 

Case 2.0<« <1. By (13), the function ¢(¢) takes values of 
opposite signs. Hence it possesses at least one zero, say t = ty. Then 
d(t)) = 0, and accordingly, by (13), 

o(t,+b) = 0. 

Hence (t+ 2b) = 0, 
and so on. Thus ¢ possesses an unending sequence of zeros at 
intervals b. (It may have more than one such sequence.) As ¢ varies 
from t, to t, +6, 4(t) varies from 0 to 0, and t+6-+-4(t) varies from 
to +6 to ty+2b. Now t+b+¢(t) is a steadily increasing function, for, 
using the variable 6 again, here {1+¢’(t)}{«+¢'(@)} = « and so 
1+-¢'(t) is always positive.* Hence as ¢ varies from t,+6 to t,+ 2b, 
¢@ passes through the same sequence of values distorted in abscissa 
and reduced in ordinate in the ratio —«:1. The behaviour of ¢ is 
therefore damped oscillatory, with quasi-period 2b, and ¢(t) > 0 as 
t > oo. Repeated application of the relation ¢’(@) -Kh’ (t)/{1+-¢4'(t)} 
shows that ¢’(¢#) > 0 as too, and hence L(t) #. The system is 


* Argument due to Titchmarsh. 





THE ENERGETICS OF NON-STEADY STATES 269 
therefore stable, the disturbances die out, and the system returns to 
the steady state. Ultimately the total energy emitted, namely, f(é), 
exceeds that generated, namely, Ht, by Ha = E(c—y)/(1+«). This 
depends on the initial conditions, as is to be expected. 





Case 3. x > 1. The argument given above shows that ¢(¢) has an 
unending sequence of zeros at intervals 6, and that the graph of ¢ in 
any interval between zeros is similar to that in the preceding interval, 
distorted in abscissa and increased (numerically) in ordinate in the 
ratio —«:1. The oscillations of ¢ therefore increase in amplitude 
without limit, and the system is unstable. Eddington has called this 
state of affairs ‘over-stability’; the system always tends to return to 
its steady state, but overshoots the mark with continually increasing 
violence. Of course in this case the physical derivation of (13) ulti- 
mately ceases to be valid, but we are now considering (13) on its own 
merits. This is the usual procedure in questions of stability. 

Case 4. « = 0. Here 7 = y = constant, and the state is steady. 

Case 5. 0 >K > —1. Put « = —k,0<k <1. Then 


d{t-+-b+o(t)} = k¢(t). (13’) 
Kither ¢ takes the value zero for some value of its argument or it 
does not. If it does, let ¢, be a zero, ¢(t,) = 0. Then ¢ possesses an 
unending sequence of zeros at intervals b, and the graph of ¢ in any 
interval between associated zeros is a distorted copy of that in the 
preceding interval, with ordinates reduced in the ratio k:1. Hence 
A(t) +0 as too and the system is stable; the disturbances are 
damped oscillatory, with quasi-period 6 (not 2b). 

If ¢ does not take the value zero, it is always of one sign. Suppose 
(i) that A(t) > 0 and let d(t’!) = m>0. Then 4(t’+b+m) = km. 
As ¢ varies from t’ to t’-+-b-+-m, ¢(t) varies from m to km, and t+-b+-4(t) 
varies from t'/+b+m to t'/+2b+m(1+k); hence, as ¢ varies from 
t'+b+m to t’+2b+m(1+k), d(t) varies from km to k?m. Repeating 
the process, we see that 4(¢) > 0 as t > 00, and the system is stable; 
the disturbances are non-oscillatory, i.e. aperiodic or “fully damped’. 


Suppose (ii) that 4(¢) < 0 and write ¢(t) = —y(t). The functional 
equation becomes b{t-+-b—p(t)} = kyb(2). (13”) 


Now «(t) cannot take the value 6. For if so b = kb, requiring b = 0 
or k = 1, both of which are excluded. Hence either always y(t) > 6 
or P(t) <b. If p(t) >, let p(t’) = b+’ (y' > 0). Then, by 
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successive application of (13”), 
Y(t —7') = k(b+7’), 
p{t’ +b(1—k)—7'(1+k)} = k*(6+7’), 
Ut + 6(2—k—k)— 9! (1+ k+14}} = 180+), 
and so on. We find eventually an argument of ys, necessarily greater 
than t’, for which ¢% is as small as we please, contradicting (¢) > b. 
Hence p(t) <b. Now let d(t’) = m <b. As t varies from ?’ to 
t’+b—m, u(t) varies from m to km, and t+b—,(t) varies from 
t’+b—m to t’+2b—m(1+k); hence, as ¢ varies from t’/+b—m to 
t’+2b—m(1-+-k), (t’) varies from km to k?m. Since m < 6b, the argu- 
ments of % for the successive intervals so constructed tend to infinity, 
and % > 0. The system is therefore stable, and the disturbances again 
aperiodic or ‘fully damped’. 

Case 6. k 1. Here (11) gives no determination of a, and the 
reduction to (13) breaks down. We return to this later. 

Case7. « << —1. Put« = —k,k > 1. The equation is (13’) above. 
If d(t) can take the value zero, a sequence of zeros at intervals b 
exists, and the usual argument shows that in each interval the 
graph of ¢ is a distorted copy of that in the preceding interval, 
increased in the ratio k:1. The system is therefore unstable (over- 
stable) and the disturbances oscillatory with quasi-period 2b. If ¢ 
never takes the value zero, it is wholly positive or wholly negative. 
If ¢ is always positive, the amplitudes of ¢ in the successive intervals 
constructed in the usual manner steadily increase; the system is 
unstable, and the disturbances non-oscillatory. If ¢ is always nega- 
tive, the usual procedure gives no result as the arguments of ¢ 
ultimately become negative. Thus the usual procedure does not pre- 
dict the ‘future’ behaviour of ¢. It is clear, however, that ¢ cannot 
tend to a finite limit as t > 00; and if it were bounded for, say, t > NV, 
then by starting with a value of ¢ sufficiently large (¢ > N’ > N) we 





could determine arguments of 4(¢) exceeding NV for which ¢ was out- 
side the supposed bounds. Hence ¢(¢) cannot be bounded as ¢t > 00, 
the system is unstable, and the disturbances non-oscillatory. 

We return to Case 6, x = —1. Equation (10) reduces to 


P{t+y—at$(t)} = d(t)+y—e. 
Choose a = y. Then 


${t+4()} = d(t)+y—e. 
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If y ~c, 4(t) cannot tend to a limit as ¢ > 00; and it can be shown 
that it cannot be bounded. The system is unstable, and as ¢(t) cannot 
take the value zero, the disturbances are non-oscillatory. 


10. We have assumed b + 0, on physical grounds. It seems worth 
while, however, to complete the mathematical discussion by a brief 
consideration of the case b = 0. 

Case 8. b = 0, —1<« <1. The equation is 

x $(t) +4{t-+-$(0} = 0. (13) 
If d(t’) = m, $(t’--m) = —«m 
d{t’ +m(1—x)} = x*m 
p{t’ +m(1—K«+x?)} = —x3m, 
and so on. Proceeding to the limit, we have 
, t’ 
ta" 
for arbitrary t’. Hence either ¢(t) = 0 or 


p(t) _ 
t+ ia constant. 


If d = 0, the state is steady. If the alternative holds, ¢’(t) is a con- 
stant, and the system has a constant luminosity L(t) = —E/«, which 


is not equal to ZH. The system therefore steadily gains or loses energy. 

Case 9.b = 0, «x > 1 or x < —1. ¢(t) oscillates with increasing 
amplitudes or tends to infinity. The system is unstable. 

Case 10. b = 0, « = 1. We have already seen from Case 1 that 
as b + 0, the amplitudes tend to zero, the periodic graph being com- 
pressed between two approaching horizontals. Thus, if « = 1, as 
b > 0, A(t) > 0 for all t. However, when actually b = 0, the equation 


becomes dit+(t)} = —d(t), (13) 
and new solutions spring into existence. If ¢(t’) = m, d(t’/++-m) = —m 
and ¢4(t’+m)—d¢(t’) = —2m, so that the slope of the chord joining 
{t’m, d(t’+-m)} to {t’, d(¢’)} is constant and equal to —2. It can be 
shown that in general the functional equation does not suffice to fix 
the ‘future’ behaviour of ¢ beyond a certain point. But an interesting 
particular case is the exact solution ¢(t) = —2t+C, which gives, 
however, ¢’'(t) = —2, L(t) = E(1—2) = —E, a case physically im- 
possible. The important point for our purpose is that (13'Y) does not 
yield periodic solutions. 
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Case 11. b = 0, «e = —1. This is equivalent to the sub-case y = c 

of Case 6. The equation is now 
P{t+-4(t)} = d(C). (13") 

Obviously ¢(¢) = constant is a solution. If ¢ is not constant, let it 
take the values A, B, (A > 0, B > 0, A # B) for values ?’, t” of t. 
Then 

g(t’+ A) = A, d(t’+2A) = A, say d(t'’+mA) = A, 

¢(t’+B)=B, ¢(t’+2B)—B, .., (t’+nB)=B. 
Now clearly A and B may be chosen so that they are not linearly 
dependent,* and then it is possible to find integers m, n so that 

(t’+-nB)—(t’+mA)| < e, d(t”+nB)—d(t'+mA)| = |A—B', 
where « is arbitrarily small. Now? if 6 = t+ (i), 
¢'(t) 
L+¢'(t) 
Hence ¢’ cannot take the value —1. But if always ¢’ < —1, then 
¢'(@) > 0, a contradiction. Hence always ¢’ > —1 and hence 
¢'(@) = 1—{1+ ¢d'(t)}}" _ < 1. Now @ passes through all large enough 
values. For otherwise |@| < constant, |4(@)| < constant, and so 
d(t)| < constant, or @ is not bounded, which is a contradiction. 
Hence for all ¢ not too small, —1 ¢'(t) <1, and so 


$'(9) 


| ts | 
$(t2)—$(t,)| = | | $'(#) dt| 
le, 


| 


Choose t, = t”-+-nB, t, = t'+mA, then 
(t,)—(t,) A—B|<e, 
a contradiction. Hence, if ¢ takes positive values, it can only take 
one positive value, and so from its continuity it is a constant. 
When ¢ takes only negative values, the argument is more difficult. 
[ am indebted to Professor Titchmarsh for the following discussion 
of this case. Put d(t) = —d(t), X(t) > 0, t>0. The equation 


aie HW} = HO. (13%) 
Case (i): #(t) 0 is a solution. Case (ii): if %(¢) > ¢ for all t, £{t—y(t)} 
is not defined, and the equation furnishes no further information 


* i.e. no integers p, gq can be found so that pA—qB = 0. 
+ The inequality now to be obtained is due to Titchmarsh. 
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about %. The system is wnstable. Case (iii): Suppose there is a value 
of t, say t,(> 0), such that 0 < #(t,) < t,. Let t, = t,—(t,). Then 
te) = Yt). If ty > W(t) = Y(t), let ts = t—P(t,). ‘Then 


ob(ts) = ob(ty) — P(t), 


and so on. We ultimately arrive at a ¢,, such that ¢,, < #(t,). Then 
the graph of ¢ lies partly above and partly below ¢(¢) = t, and hence 
there exists a ft, such that #(t)) = t), whence by (13%) %(0) = ty. 
Since 7% is single-valued, there is only one ¢t). The graph therefore 
passes from (0,t)) to (f),f)) and for t > tf, lies below y(t) = t. It is 
clear that any value taken by (i) anywhere is also taken in (0, t,). 
Hence (tf) is bounded. By considering the intersections of the 
graph of %(¢t) with the lines 4 = t/n, Titchmarsh now shows that 
these intersections only occur where % = ft), and hence that, for 
nly St < (n+1)ty, n(n+1)t, < A(t) < (1+n)t,. In particular, 
s(t) > ty as too. It can then be shown that as ¢, varies from nt, to 
(n+-1)to, t varies from (n—1)t) to nty, tz from (n—2)ty to (n—1)to, 
and so on. Moreover, since dt,/dt, = 1—y'(t), and since y#'(t), by a 
modification of an argument previously used, can be shown to be 
everywhere < 1, fg, ty,... increase through the intervals mentioned. 
Hence ¢(t) passes through the same values in each of the intervals 


{nty, (n+ 1)tp}, and since b(t) > ty as t > 00, #(t) = ty always. Hence 
L(t) = E, and the state is steady. 


11. We have now shown that the functional equation (13) generates 
of itself, without detailed assumptions about the nature of the system 
under discussion, all the different dynamically possible types of dis- 
turbance, and that the only case in which strictly periodic solutions 
occur isk = 1,6 #0. That is, if the supply of energy is at a constant 
rate H, periodic fluctuations in the emission can only occur, by 


(7), when (22) = 1 
o-0, 


ea _ 


9 
dQ Sa 


12. Application to Cepheid variation 

It follows that a star can only* be a Cepheid when very special 
conditions are satisfied. In general, for a star which exists, « will not 
happen to be unity, and disturbances will rapidly damp out. For an 


* Assuming that the rate of generation of energy is unaffected by dis- 
turbances and that T is a single-valued function of 2. 
3695.4 T 
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ordinary star, then, —1 < « < 1, or 


dT 
—1l< —H|—~ i 24 
a (ia), 2 = . ( 


the central term of the inequality is a sort of damping parameter 
(dimensionless) which determines the rapidity with which the oscilla- 
tions die away. When the damping parameter nears the value —1, 
the star begins to be unstable. When it nears the value +1, the 
oscillations tend to become of constant amplitude and strictly 
periodic, though asymmetric in form. In the limit when the damping 
parameter is unity, undamped periodic oscillations are possible. If 
we call this the Cepheid state, it follows that the Cepheid state can 
continue indefinitely. It has long been a puzzle, (a) why all stars are 
not Cepheids, (6) why Cepheid oscillations maintain themselves 
apparently undamped. Our theory answers both questions. Only in 
particular circumstances will a star satisfy condition (23), so that the 
Cepheid state is not in general possible. But if condition (23) is 
satisfied, strictly periodic undamped oscillations are compatible with 
a constant rate of generation of energy. Further, the theory predicts 
that the fluctuations in the energy output at the surface will not be 
of simple harmonic form, and that maxima will be more pronounced 
than minima, in agreement with the observed characteristics of 
Cepheid light-curves. 

The quantitative relations (19) and (20), or (19’) and (20’), together 
with (21) and (22) offer an ample field for further tests of the theory. 
Whether they are found to be obeyed or not, the analysis suggests 
that attention could with advantage be fixed on the relations between 
the asymmetry of the light-curve (on an absolute energy-scale) and 
the amplitudes of the oscillations measured from the mean ordinate. 
Cepheid light-curves have usually been discussed with reference to 
the median ordinate of the magnitudes. Further, the asymmetry 
should be measured by our quantity 5,* the interval between the 
attainment of mean luminosity on the upward branch and that on 
the downward branch, and not in terms of the interval from maxi- 
mum to minimum. ‘Steepness-period’ relations have been found by 
Ludendorfft and Vernon Robinson,{ but they are not easily related 


* Better, by 5/P. + Berlin Sitz. v (1929). 


t Harvard Obs. Bull. (1930), 872 and 880. Cf. K. Lundmark, Lund Obs. 
Circular, No. 7 (1932), 150. 
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to our predicted relations between mean amplitude and asymmetry 
of energy-curve. 

If relations (17)—(20) are not found to be obeyed by actual Cepheids, 
the conclusion would be that the pulsations affect the rate of genera- 
tion. 

Our theory shows immediately why the Cepheids form a definite 
sequence in the Hertzsprung-Russell diagram or any equivalent plot. 
A star is defined by two independent variables, luminosity ZL and 
mass M, assuming a constant mean chemical composition. Thus our 
damping parameter « will be a function of Z and M, x = x(L,M). 
The condition «(L,M) = 1 defines a locus in the plane of Z and M, 
which should be the Cepheid sequence. All the remaining tenanted 
portions of the plane should correspond to —1 < «x < 1. This shows 
again that Cepheids should be relatively infrequent amongst the 
stars. Further, the parameter b = (y+-cx)/(1-+-«) will also be a func- 
tion of L and M. When « = 1, P = 2b, and between the relations 
«(L, M) = 1, b(L, M) = $P, we can eliminate M, obtaining a relation 
F(L, P) = 0. This is the well-known period-luminosity relation. 

It is quite unknown whether the stars evolve or not. If they do, 
a Cepheid stage would only be reached when « secularly increased 
up to unity; it would be passed when « began to decrease again. It 
is clear that though a star would be unstable for a value of « not 
lying inside —1 < «x < 1, it is impossible to have an actually unstable 
disturbance occurring except for values of « infinitesimally greater 
than unity. For a star could not reach a value of « not lying inside 
—1< x« <1 without passing through unstable states, which is im- 
possible. 

A further test of the theory would be a direct physical calculation 
of T(t) or 7. For we have seen that 7’ = b = }P, so that the observed 
period of a Cepheid should be twice the mean value of the calculated 
time-lag. We have seen that b = T = }(7,+7,), so that the mean 
lag T is not necessarily equal to the ‘steady-state lag’ 7;, though it 
must be admitted that the ‘steady-state lag’ 7, is a somewhat con- 
ventional quantity, incapable of direct calculation from its apparent 
definition. In connexion with this it may be noted that, from (2) and 
(9), the mean energy imprisoned in the system is 


Q = Q,— Ea = Q)—$E(c—y) = Q)—$H(T,—T)). 


Now 6 and Q must be independent of the choice of origin of time, 
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whilst 7, and Q, depend solely on the choice of the origin of time. 
Hence 4$(7,+7,) and Q,—4}#(7Z,—%) are independent of the choice 
of time-origin. It follows that the conventional quantity J, asso- 
ciated with a given periodic disturbed motion depends on the choice 
of time-origin. If we choose the time-origin so that Q, is equal to 
the steady-state energy Q,, then 7, = 7, = 7. This is physically 
satisfactory. 

13. Application to ‘Griffiths’s Heat Engine’ 

In this apparatus the bend of a U-tube is sealed with mercury, 
and one arm of the U sealed to a flask containing air, the other arm 
being open. When the flask is steadily heated, the mercury begins 
to oscillate, the oscillations being reported to be of very considerable 


amplitude. The assumptions of our general theory seem to be ful- 


filled, probably by adjustment of the rate of supply #. It will be 
noted from (7) that increase of # increases x numerically, and so 
shifts the system in the direction of reduced stability. This is physi- 
cally obvious, for clearly sufficient increase of H will make the engine 
unstable. It would be highly interesting to test experimentally 
whether the oscillations are of the asymmetric non-harmonic type 
predicted by our theory. This seems probable, from Griffiths’s own 
discussion,* for he points out that the average air-pressure in the 
flask must be greater during the descent of the mercury (on the flask- 
side of the U) than during the ascent. The energy supplied to the 
flask is dissipated by radiation and conduction from the tube. 

14. Functional equations of the type (13) make their appearance 
in other physical problems, in particular in the theory of the syn- 
chronization of watches by two relatively stationary observers by 
means of light-signals. It would render the present paper too long to 
discuss this problem here. The broad point of physical interest in our 
analysis is the proof that periodic fluctuations of a system are, under 
certain circumstances, compatible with a constant supply of energy, 
and that such fluctuations are necessarily not simple harmonic. 
Note added in Proof. 

Mr. J. Hodgkinson has kindly pointed out that, by use of the 
method of finite differences,t formal solutions of the equation 

Kb )-+4{t-+b-+4(t)} = 0 


* Loc. cit. 51. 
+ G. Boole, Calculus of Finite Differences (2nd edition, 1872), Ch. XV. 


. 
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may be obtained in the following forms: 

(a) ck >0; f(t) = C(x)x*sin rx, 
where t = A(x)+ba—C(x) — ff 

I1+k 

(b)x <0, n = —k (k $1); 

P(t) = C(x)ke, 
where t= A(x)+bx+ Oa); 
(c) «= —1; o(t) = C(x), 
t = A(x)+br+2C(z); 
where in each case C(x) and A(x) are arbitrary functions of period 1. 
Many of the properties of ¢(¢) obtained above can be derived from 
these explicit solutions. 

The paper was, however, concerned with functions ¢(¢) which were 
single-valued functions of ¢, and Mr. Hodgkinson’s solution does not 
yield a single-valued function of ¢ unless some restrictions are placed 
on C(x) and A(x). Consider, for example, the particular case x = 1, 
for which the solution is 

f(t) = C(x)sin 7x 
t = A(x)+ba—$C(x)sin ax. 

A(t) is periodic of period 2b, and x increases by 2 as ¢ increases by 2. 
But if it ever occurs that dt/dx becomes negative, the graph of ¢(t) 
will turn backwards and ¢ will cease to be one-valued; for the graph 
of 4(¢) can be obtained by plotting C(x)sin 7x against bx and shifting 
each ordinate forward a distance A(x)—}C(x)sinzx. In particular, 
as ¢t varies from ¢t to t’ = t+b+¢(t), x increases by 1, and accordingly, 
in order that ¢(¢t) shall be single-valued, it is necessary that 
i(x+1)—t(x), the change of ¢, shall be positive, i.e. that b+ C(x)sin 7x 
shall be positive, i.e. that |C(x)sinax| < b for all x. The inequality 
—b < d(t) < b then follows. If C(x) is capable of such a value C(x’) 
that C(x’)sina7x’ = —b, then t(x’+1)=t(2’), but d(t) = —), 
d(t’') = +6, and ¢ is not single-valued. 


where 
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MEAN-VALUE THEOREMS AND THE RIEMANN 
ZETA-FUNCTION 
By A. E. INGHAM (Cambridge) 


[Received 12 September 1933] 


1. THE Riemann zeta-function {(s) = {(o+ti) satisfies a mean-value 
theorem of the form Lr 
. 1 y } ° 9} . 
lim 7| \f(o+tz) |?* dt = C(c) . (1) 
T—>2 
in the following cases: 
(I) when & is real and o > 1; 


: ash ae 1 1 
(Il) when & is a positive integer and o > max -, 1——]; 
£ > 
- v 


(III) when & is a positive integer and. o > }, if the unproved 
‘Lindeléf hypothesis’ concerning ¢(s) is true. 

Of these results, (I) follows (in virtue of the expansion of (*(s) 
given below) from a classical theorem of Hadamard, (II) is due to 
Schnee in the case k = 1 and to Hardy and Littlewood in the case 
k = 2 and may be extended to k = 3, 4,... by the methods of the 
last-named authors (or by a method developed by Carlson for the 
case k = 2), and (III) is due to Hardy and Littlewood.* 

The question of the existence of the mean-value (1) for positive 
non-integral values of k and for a range of values of o < 1 has been 
considered by Titchmarsh.t He has extended completely the result 
(III) which depends on the Lindeléf hypothesis, and, as an example 
of what may be obtained from his methods without unproved hypo- 
thesis, he has shown that the mean-value exists for k = } and o > }. 
In this paper we shall extend the result (Il) completely to positive 
non-integral values of k; thus, in particular, we show that the mean- 
value (1) exists for k = } and allo > 3. 

In all these theorems “en 

C,(c) = S din) 


n2o 
n=1 


* For details and references see Landau (7), (ii) 776-8, 806-19, 905-6, and 
Titchmarsh (9), Ch. I (1.23), II, VI. (The references in heavy type are to 
the list at the end of this paper.) 

+ Titchmarsh (8). 
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where the coefficients d,(m) are defined by the following expansion, 


which is absolutely convergent for o > 1: 


-T¢-3) = [+e +) - 2 


The series an C,(c) is convergent for all o > 3 1 ne ; 
d,(n) = O(n*), 
for every fixed positive «, when n > 00.* 
The proof of our extension of (II) is based on (I) for k > 0, on the 
special case k = 2 of (IL) (or rather the weaker theorem that 


T 
7 | Ke+eii¢ at = 00) 


as T' > o, for every fixed o > 4), and on the Bohr-Landau theorem 
that, if 5 is fixed and positive, the number of zeros of {(s) in the 
region o > $+8, —T <t< T, is o(T) as T +o. From these and 
simple classical properties of {(s) we deduce our theorem by the 
application of certain theorems of general function theory. These 
theorems concern mean-values along lines, and are extensions or 
adaptations of theorems considered previously by Carlson, by Hardy, 
Ingham, and Pélya, by Gabriel, and by Besicovitch.t We state and 
prove these general theorems in § 2 and 3, and apply them to {(s) 
in §§ 4 and 5. 

2. Theorem A. Suppose that f(s) = f(¢+ti) is regular in the half- 
strip S defined by 0, <0 < ag, t >to, and satisfies a ‘Phragmén- 
Lindeléf condition’ 


f(s) = Oe"), 0<h< 


Og—9, 


uniformly in S. Suppose that (when T > o) 


T 7 
lim > ( [ \flo,+ti) " a <a 


lim re | o+tt) |N% a)” « A, 


* This is classical when k is a positive integer, and the general case follows 
from this since |d,(n)| < dn) where J is an integer greater than |k|. [Cf. 
Titchmarsh (8).] 

+ Bohr and Landau (2). There is a more precise theorem due to Carlson, 
but this and its further refinements are not required. [See Titchmarsh (9), 3.6.] 

t Carlson (4); Hardy, Ingham, and Pélya (6); Gabriel (5); Besicovitch (1). 
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where x1, %, 44,49, A,,A, > 0. Let 


3, > 0, 3, > 0, 3,+8, = 1, 


and let o, x, a, A be defined by 
o = 3,0,+%,0,, a= B,a,+0,a,, a= ,a,+80,a, A= A? 42%, 


T x 
Then lim re( f(o+ti) |" at <A, (3) 
fo 


and this‘holds uniformly in 3,0, in the sense that, given any positive 
e, we can find a T, independent of 3,,3,, such that the expression 
following the sign lim in (3) is less than (1+«)A for all T > T, and 
all admissible values of 3, 35. 

This theorem is an analogue for mean-values of Theorem 2 of 
Gabriel’s paper 5. It is related to Theorem 8* of Hardy, Ingham, and 
Pélya’s paper 6 as Gabriel’s theorem is related to Theorem 7 * of that 
paper, and (apart from minor complications arising from the use of 
unilateral means) may be deduced from Gabriel’s theorem in much 
the same way as Theorem 8 was deduced from Theorem 7. We begin 
by proving two lemmas which embody the main ideas of the proof. 

LemMaA 1. Suppose that f(s) 1s regular in the rectangle R defined by 
OS 0S 0g, 4, <t < T; that 

T Oy T Xs 
( | | f(o,-+d2) [No a < A,, ( | | f(ea+-#i) | i) i.e oe 
te to 
where x1, %, Ay, Ay > 0; and that 
|f(s)| <n min(A,, Ay) (4) 


on the horizontal sides t = tj and t= T of R. Let 31, by, 0, a, A be 
as in Theorem A. Then 


m 
( [ | f(a ta) |! at) < HA, 
F 


where H is a function of 04, a2, x, %, n, which tends to 1 when n > +0 
(1, Go, 4, X remaining fixed). 

Denote the left-hand side of (5) by /(3,,3,) or J(%). We first con- 
sider pairs {3,, #,} of the form {r,/2”, r,/2"}, where r,, 72, n are integers; 
such a pair will be said to be of order n if one of the fractions 7,/2", 
r,/2” (and therefore also the other) is in its lowest terms. Let n > 1, 


* Or the corresponding theorems of Carlson’s note (4) (398). 
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and suppose we have proved that 
I(t) < H,_,A 
for all pairs {3,,%,} of order not exceeding n—1. Let 
{91,8} = {r1/2",79/2"} 
be any pair of order n, so that 7,, 7, are odd. Then the pairs 
(9), OS} = {(7, +1)/2", (r2—1)/2"}, (8, 83} = {(r1—1)/2", (r+ 1)/2"} 
are of order not higher than n—1, so that, by our supposition, 

(8) <H,44’, 1) < 0,44", (6) 
where o’, «’, A’ relate to the pair {3{, 3}, and o”, a”, A” to {H, BH}. 
Now, since 3, = 4% +4 and 3, = 48;+ 48), we have 

o = }o’+}o’, a == $0’+ 40”, A= A'tA"); 
also c—o’ = o”—o = (o,—o,)/2". Hence, by a theorem of Gabriel,* 
applied to the rectangle bounded by the lines of abscissae o’, o”, and 
ordinates t), 7’, we have 


J ineyn asl) *{ f isteye asi) 
z 


K’ 


1(8) < 


where K’, K” are the two halves of the boundary of the rectangle 
which lie to the left and right respectively of the line of abscissa co. 
Now, by (6) and (4), the first factor on the right of (7) is less than 


((H,-1 A’)Y%+2(nA’)"*(ag—03)/2"}* 
(HES + (oy —oy)2-Dqhlapiora” 
(HY, + (oq—0,) 272-0 Datty glint A 
(Hi ly rs Uu Aply)ivA’s, 
where 


y = max(ay, a9), U = max (o,—¢,)””, u = max 2-7; 
X= My ,Qy X= My yg 


we have used here the inequality 2+y < (a«+y) («# > 0, y > 0, 


* Gabriel (5), Theorem 1. The need for a result of the type of Lemma 1 
arises merely from the use of unilateral means. If we were working with means 
over (—T’, T) in Theorem A [which would, of course, suffice for the application 
to {(s)], we could simply quote Gabriel’s Theorem 2. 

The characteristic feature of Gabriel’s Theorem 2 is convexity in two 
variables (o and « in our notation). Other theorems embodying this feature 
were proved later by Carleman (3), and there is an earlier remark of Carlson 
asserting without proof a theorem of a similar nature [Carlson (4), near the 
bottom of p. 397, ‘pour « > 21’expression v/x est une fonction convexe des deux 
variables o et 1/x’]. The results given by Carlson and Carleman are insufficiently 
precise for our purpose. 








282 A. E. INGHAM 
A> 1), with A = y/«’. A similar result holds for the second factor 
on the right of (7). Hence 

I(8) < (HY, + U0 un )\vAt A"? = H, A, 
where Hi’y = Hy, + U0u-1y). 
Now for order 0, i.e. for the pairs {1, 0} and {0, 1}, we obviously have 
I(d) < A. It follows by induction that, for all pairs {3,, 3,} of order 
not exceeding n, /(#) < H, A, where 


n 
Vy — Tyr-lyly — 
LY = ow nly < 


since 0 < uw < 1. Hence (5) holds, with 


U ie 
H = ( -— it) : 
l—u 


for all {3,, 9,} of the special form considered, and therefore, by con- 
tinuity, for all {3,,3,}. Since H depends only on oj, a9, a, a), 7, and 
tends to 1 when 7 > +0, the lemma is proved. 


Lemma 2. Let E(t) = E,(t) = e-“, where kk > 0. Then 
r 
he vr | p(t) dt = im ~ | E(t—w)d(t) dt (8) 
for any iil (integrable) 4(t) and any real c. Further, if 
T > max(ty, 7 and w = T'+log T’,, then 


ak at < (14 Fe) er FF Gl 


Denote a left- and right-hand sides P (8 3) by A, and A, respec- 
tively, so that 0 < A,, A, < +0 (with obvious conventions about 
+00). 

We first prove that A, > A,. We may suppose that A, < +00 
since otherwise there is nothing to prove. Then we can choose num- 
bers C > 0 and r > max(fy, 0), so that 

r 
o(T = J #0 t) dt < CT° (T > 7). 


Since 0 < E(t) < 1 for all t, and d(t) > 0 for t > ty, we have 
oa) i 
| Bi—w)$(0) dt < | at + f B(t—w)p(t) ae (TF > &2. 


to to i 








iy 
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Integrating by parts, and noting that —H#’(t) > 0 for all ¢t, and 
0 < O(t) < Ct for t > 7, we see that, if 7’ > 7, the last integral is 


ao 


— E(T—w)®(T)— | h (t—w)@(t) dt < — | E'(t—w)Ct dt 
T > 


2 ~ e 

| E’(u)C(u+w) du = — | BwyCer(*+1) du 
Tw T-w “s 

Hence, if 7 >w>1, 


Ls | 2e—aype at < (7 s) Fe | 40 t)dt — | B’(u)C(u+1)" du. 


Ww 


« 


to 
Now let 7 and w tend couttuniesai to +00 in such a way that 
T/w—>1 and T’—w > +0); we infer at once that A, < A,. 
On the other hand, since H(t) is positive and decreasing, and 
d(t) > 0, we have, for w > 0 and T' > ty, 


oO i 
5 [ E(t—w)d(t) dt > (=) Bro) re | d(t) dt 
t, ‘4 


Making 7 and w tend continuously to +00 so that 7'/w—1 and 
T—w—> —oo, we deduce that A, > A,. And taking, in particular, 
- T+log 7’, we obtain the last part of the lemma.* 


Proor oF THEOREM A. Choose k so that 


A<bh< ey 
Og—90, 


and let 
e(s) = ee Kt! (8) = e(s—wit) f(s) (w > 0). 


In the half-strip S, 


le(s)| = e-eMcos k{o—Ho,+0%)) 
1 
2 


and 1 > cosk{ao—}(o,+-04)} > cos $k(o.—o,) > 0, 

* Lemma 2 is modelled on Hardy, Ingham, and Pélya (6), Lemmas B and C 
(557-9). The latter might be used instead, but Lemma 2 seems a little simpler 
because it avoids the additional limiting process m-— oo, and airc more 
natural because a function of the type of Z(t) must in any case enter in con- 
nexion with the Phragmén-Lindeléf condition. 

In stating Lemma 2 we have confined ourselves to what is actually required 
for the application. But we may add that, if A, = <A, is finite, then the lower 
limits A, and A, are also equal (as a glance at the proof will show). We thus 
have a (trivial) Abelian-Tauberian relationship. 
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since 0 < 4k(c,—0,) < $a. Hence we can choose a positive number 
p (depending on oj, o9, a4, %, k) such that 
ee < a cos k{a—4(o,+02)} < ek 
for 0, < o < o, and all a between «, and a,; and then we shall have 
E(t+p) < |e(s)|"* < E(t—p), 
where EL E(t) has the same meaning as in Lemma 2. Hence 
E(t—w+p)|f(s)|¥* < |f,,(s)|"* < E(t—w—p)| f(s) |! (9) 
for all s in S and all « between a, and ag. 
Now by hypothesis 
. 
1 
, . | fo\|1/a < 1/Qy 
bw Tale, | f(a, + ta) |}! dt < Al™, 
to 
Hence, by Lemma 2 (with w+-p written in place of w), 
RD 
: ] * : . : 
lim ie [ E(t—w—p)|f(o,+t) |" dt < Al, 
wre (wp) J 


to 
This implies, in virtue of the second inequality (9), 


ao-rn 


c— Ulfis 4 
lim =i J | Fu, +-e2) |" at <A, 
to 


[since (w+p)/p > 1]. Thus, if any positive number 7 is given, we 
can find an w, = w,(7)* such that 


( j froloy-+ti) at)” <fivics eal 
t 


and similarly an w, = w.(7) such that 


( [ | fos tt) |N% at ss < (1+ n)A,w® 
to 
Now, if M is the maximum of | f(s)| on the segment t = fy, 0, < o < og, 
we have, on this segment, | f,,(s)| < |f(s)| < M and therefore cer- 


tainly f..(8)| <7 min{( 1+ )A w, (1+7)A,w® 


for w > w’ = w'(n). And for any fixed w the same inequality must 
hold on the segment t = 7’, o, < ¢ < og, for all sufficiently large 
T'- for 
f(s) = Olexp{—e e-*cos hk(o,—0,)+e"}] = 0(1), 
* We do not indicate explicitly the dependence of w,, ws,... on ‘constants of 
the problem’ such as f, ty, 01, G2 4, XM, Aa, Ay, Ay. 
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uniformly in S, ast > 00, since k > h > Oand e-*“cos $k(a,—<a,) > 0. 
Hence, taking & > max(w,, w, w’) = wy, applying Lemma 1 to f,,(s) 
and the rectangle o, < o < og, tj < t < T, and making 7’ > 0, we 
obtain 

- a 

(J \fa(o-+ti) a) <H(l+n)40" — (w > @), 
to 
where o, a, a, A are as in the enunciation of the theorem, and 
H = H(y)>1 when »—> +0. Using now the first inequality (9), 
changing w to w+ p and dividing through by w*, we deduce that, if 


w > wy (> &—p); 


( 
w"\ J 
to 


It now follows from the last part of Lemma 2 (with c = a/«) that, 
if T > max(ty, wp», 1) and w = T'+log T (> w,), then 


we 
- a 


E(t—e)| flo-+ti) at)" < H(1+7)4(1+2) 
Ww 


1 r a p a log T a = 
| | | f(a+ti) |" a < H(+9)4(1+2) (1+ 2 ) er. 
[a ; W 7 
Now, if y = max(a4, a), 9 = max(@,, 4), the product of the last three 
factors on the right is 

p\" log T\9 p-* 

PY Ty y 1+ 

< (145) (14%, et <14+7 

for T > T’ = T’(n). Hence, if 7’ > max(to, wo, 1, T’) = T% = T(n), 
1 


l - 
7a | 
to 


Now, given « > 0, we can choose 7 = (e) so that H(1+7)? < 1+e, 


|flo+ti) at)" < H(1-+n)*A. 


and then the corresponding 7, = 7)(y) = 7 (e) fulfils our require- 
ments. 

3. Theorem B. Suppose that f(s) is regular in the half-strip S 
defined by 0, < o < ag, t > ty, and satisfies the following conditions: 


(i) for a given q > 0 the limit 


a 
lim 7 | Iflo+-ti)|@ dt = M,(0) 
Tox 

to 


exists for an infinity of values of o having a limit point inside the 


interval a, <a < G9; 
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(ii) for some Q > q 
T 
7 [ Iflo-+ti)|@ dt = O(1), 
te 
uniformly for o, < o < ag, when T > ~; 

(iii) N(T) = o(T) 
as T > 0, where N(T) denotes the number of zeros of f(s) in the rect- 
angle o, <a <on,t, <t<T. 

Then the limit M,(c) exists for all o in the interval ao, < o < og, and 
uniformly in any interior interval o, < 0, < o < a3 < a9. 

Further M,(c) is the value for s = o of a function M,(s) = M,(o+t) 
regular in the strip 0, <0 < oy. 

When q is an even integer, no restriction on N(7’) is needed, and (ii) 
may be replaced by the corresponding condition with Q = q (a wider 
hypothesis in virtue of Hélder’s inequality); this (or rather the ana- 
logous result for bilateral means) constitutes Theorem 14 of Hardy, 
Ingham, and Polya’s paper 6. The result was extended to general 
q by Besicovitch, subject to the hypothesis (iii) and to a different 
type of condition in place of (ii).* It is the substitution of (ii) for 
the corresponding condition of Besicovitch that is important for our 
application to {(s). The condition (iii) will be satisfied in the applica- 
tion, in virtue of the theorem of Bohr and Landau quoted in § 1. 
For convenience we give the proof of Theorem B in full, though it 

* Besicovitch (1), 375. In our notation Besicovitch’s condition takes the 
form 


T+1 
(ii)’ f |f(o+ti)|¢dt < C (o, fe <o, T>k). 
; 


But Besicovitch considers, in place of our functions f%, functions which may 
have singular points of any kind in S, provided that their distribution is 
subject to a condition of the type (iii) and that (ii)’ is not violated. Thus our 
theorem neither includes nor is included by his. For the special class of func- 
tions considered here, however, the substitution of (ii) for (ii)’ represents a 
generalization, since (ii)’ in S implies (ii) (and indeed the boundedness of /) 
in any interior half-strip Ss; of the form o,+6 < o < o,—8, t>%+8. In 
fact, if s* = o*+t*i is any point of Ss, we have, since | |? is subharmonic, 

r) T o*+8 #46 

78" | f(s*)|¢ < fr dr f |f(s*+re%) |" d0 < f da f |f(a+ti)|* dt < 28C, 

0 7 o*—5 *—§ 
if 8 < 4 (as we may suppose without loss of generality). This type of argument 
(which has been employed by various writers) may be used to simplify the 
proofs, and generalize the results, of Theorems 2, 10, and 11 of Hardy, Ingham, 
and Pélya (6). 
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is only in connexion with the condition (ii) that there is any point 
of novelty. 

Let the distinct ordinates of zeros of f(s) in the interior of S be 
t,, ty,..., Where fj << t, << t, <.... Take 1 > 0, and denote by L(7’) 
the set of points ¢ of the interval (f), 7’) lying outside all the intervals 
(t,, —l, t, +1) (n = 0,1, 2,...), and by 1(7’) the remainder of the inter- 
val (t), 7’). If 1(7') denotes also the measure of the set 1(7'), we have, 


by (iii), UT) = 0(T) (10) 


as T’'+> oo. Using Hoélder’s inequality, (ii), and (10), we now obtain, 
writing f = f(o+t), 


l l alQ l \(Q-M/gQ 
nae | f\@ 4. | f!Q peo QIQ-a) 
> [ if w<(p f if a) (7 i: i) 


UT) UT) UT) 
» 


é (7 | f/2 ee 0(1) 0(1) = 0(1) 


0 


uniformly for o, < o < oy. Hence 


Lir) 


Y 
7 | Mo+eyed — 7 [ Wo+tiie do, (11) 
ty 


uniformly for o, < o < og, when T' > ©. 

Now the lines ¢ = t,, divide the interior of S into domains in which 
f(s) is regular and has no zeros. Define a function g(s) in this set of 
domains by taking a definite branch of f!%(s) in each domain, and 
let g(s) be the conjugate function (so that g(s) = g(s) and is regular 
in the images of the above domains in the real axis). Let 


G(s, T) = Gis, T) = z | g(s+ui)g(s—ut) du. 
L(T) 


For fixed 7’ > max(ty, 0), G(s, 7’) is evidently regular in the domain 
D, defined by 0, < o < az, |t| <1. Also 


, i 
|G(s,T)| < (7 | lewd} (7 | ate—ui da 
L(T) L(T) 


T+1 
<p | Metre = 0(*-) = 00), 
to 
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uniformly in D,, as T' + 00, by (ii) and Hélder’s inequality. Further 
(if o is real) 

ir 
ae 
L(T) L(T) 


G(o, T) = LA lg(o+ut) |? du = \f(o+ui)\¢du, (12) 


- 
so that, by (i) and (11), G(s, 7’) tends to a limit, as 7’ > oo, for an 
infinity of (real) s having a limit point in D,. It follows, by Vitali’s 
theorem*, that G(s, 7') tends to a limit G(s) = G,(s) at all points s of 
D,, and uniformly in any closed set in D,; and G(s) is regular in D,. 
Hence, by (11) and (12), 


a 


i 
-_ 


, 


\f(o+ti) q dt —> G(c) 


for o, < o < o,, and uniformly in any interior interval. The last 
relation shows that the values of G(s) on the real axis are independent 
of 1. Since / may be as large as we please, it follows that there is 
a G(s) regular in the whole strip o, < o < o, taking the values in 
question on the real axis. This completes the proof. 


4. We can now state and prove our main theorem concerning {(s). 
The symbol C,,(c) has been defined in § 1. 


Theorem. Let d be a fixed positive number such that 
» 


= [ L(o-+ti)|\ dt = O(1) (13) 
J 


] 


as T' +>, for every (fixed) o > 3. Let o, and q be two fixed numbers 
satisfying the relations 
(1—o9)q < 2A. 


Then 7 \C(a-+t2) |¢ dt > C,,(o), (14) 
» 
uniformly for o > ay, when T > ~. 

We first apply Theorem A with f(s) = €(s); tj = 1; o, = $+, 
o, = 1+8, where 0 < 26 < o,—}; a, = a, = 1/A; ag =a, = 5. All 
the hypotheses are satisfied; the two upper limits involved are finite, 

* See, e.g., Titchmarsh (10), 168 (5.21). The theorem is there stated for 


a sequence of functions, but the proof given applies also to a set of functions 
depending on a continuous parameter. 
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the first in virtue of (13) and the second because {(s) is bounded on 
the line o = 1+8. Hence 

r a 
(3 [coe a) = 011) (15) 
1 
uniformly in #,, #,, where #,, #,, 0, x are as in Theorem A. It follows 
from Hélder’s inequality (since o, < og—8 < o,) that (15) will hold 
uniformly for og—8 < o < oy, if « is replaced by any fixed number 
greater than &, the maximum of the values of « corresponding to this 
range of o. Now, since #, = (o,—a)/(a,—o,) and #, < 1, we have 


_1—o9t+ B 


+A = - 1) — +6, 


a< o2—(o,—8S) 
O2—9; 
Since (1—o,)/($A) < 1/g, we can choose 8 = 8(A,o9,q) so that 
¥ < 1/(¢+8), and we shall then have 
° 
1 
(7 | 


oy 


1/(q+8) 
IK(o-+ ti) |a* a) = O(1) (16) 


uniformly for o,—8 < o < 1+. 
We now apply Theorem B to ¢(s) and the half-strip 
o—5 < o < 1-44, ey 
We know [see § 1, (I)] that (14) holds for 1 < o < 1+86; thus the 
condition (i) of Theorem B is satisfied. Condition (ii) is fulfilled in 
virtue of (16), and condition (iii) by the Bohr-Landau theorem (since 
co )—d > 4). Hence, when 7' +> 00, the left-hand side of (14) tends, 
uniformly for og < o < 14-36, toa limit M,(c), where M,(s) is a regu- 
lar function of s in the strip o,—é < o < 1+8. Now, when s = oa, 
1<o< 1448, M,(s) coincides with the function 
 di,(n) 
C,,(8) — bi <2 ? (17) 
n=1 

which is also regular for o,—85 < o < 1+8, the Dirichlet’s series being 
(absolutely) convergent in this strip*. Hence M%,(s) = C,,(s) through- 


* This is obvious from the fact that d,, (n) = O(n‘), but may be inferred 
independently. For if o, is the abscissa of convergence of (17), we must cer- 
tainly have o, < 1, since the expansion > d,,(n)n~* defining d,,(n) is abso- 
lutely convergent for ¢ > 1, and so a fortiori is (17). And, if ox—8< o, < 1, 
then, since dj, (n) > 0, C,,(s) would have a singular point at s = o, (by a 
well-known theorem of Landau), and this is incompatible with the existence 
of the function M, (8). 

3695.4 U 
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out the strip, and in particular for s = o, og < o < 1+38. We have 
thus proved that (14) holds uniformly for og < o < 1448. But it 
also holds uniformly for o > 1+ 46 (as a glance at the usual proof 


for o > 1 shows). This completes the proof of the theorem.* 


5. In the theorem just proved we may take A = 4, in virtue of 
(2). The resulting proposition is equivalent, except for the difference 
of notation and for the additional assertion about uniformity, to that 
enunciated in § 1, namely, (II) with the restriction to integral values 


of k removed. 

We may note also that, if the Lindeléf hypothesis is true, then, by 
(III) of § 1, we may take A to be an arbitrarily large even integer, 
and therefore o, arbitrarily near to } and q arbitrarily large. We 
thus obtain an alternative proof of Titchmarsh’s extension of (III) 
to non-integral k. For this we do not, indeed, require the full force 
of our theorem, for it is evident that a convexity theorem in the 
two variables o and « is not needed in this case. 


* The second half of the argument embodies the proposition that, if the 
left-hand side of (14) is bounded for a given q (and an appropriate range of 
a), then (14) is true for any smaller (positive) g. A similar result has been 
obtained independently by Bohr and Jessen in a forthcoming paper. These 
authors apply their methods to some further problems where the argument 
of the present paper would break down owing to the absence of any analogue 
of the Bohr-Landau theorem. 
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THE WEIERSTRASS E-FUNCTION IN DIF- . 


FERENTIAL METRIC GEOMETRY 
By J. H. C. WHITEHEAD (Ozford) 
[Received 18 October 1933] 


1. Iy this note it is shown how the function which is known in the 
calculus of variations as the Weierstrass EH-function* is related to 
the idea of convexity dealt with in two previous notes by the present 
author.t 
In these notes we proved the existence of ovaloid hypersurfaces in 
a space of paths, a hypersurface being described as ovaloid if it 
bounds a simple and completely convex region (1). Here we show 
that the regularity condition, H > 0, is necessary and sufficient for 
small spheres, defined in terms of the fundamental integral, to be 
ovaloid. 
The condition EZ > 0 is seen to be a consequencet of 
g = |9i;| A 9, F(x,dx) > 0 (dx + 0), (1.1) 
1 @F? 
where Jij(x,dx) = 3 ddd 
The conditions (1.1) also imply 
9i;(v, dx) Sada) > O (dx + 0, 5x # 0). 


2. Let the fundamental differential invariant, 


F(z, dx) = F(z"....,a°, d2",..., dz"), 


be of class 2, homogeneous of the first degree in dx, and let it satisfy 
the conditions 


g ~ 9, F(x,dx) > 0 
for all values of 2 in a given region and all non-zero vectors dx. The 
Weierstrass function is defined as 

E(x, p, é) - F(x, £)—F,.€, 

* For an account of this function see, for instance, J. Hadamard, Legons 
sur le calcul des variations (Paris, 1910), Livre ITT, ch. II, or O. Bolza, Vorle- 
sungen tiber Variationsrechnung (Leipzig, 1909). Also we refer particularly to 
a paper by G. A. Bliss, Z'rans. American Math. Soc. 15 (1914), 369-78. 

+ Quart. J. of Math. (Oxford), 3 (1932), 33-42, and 4 (1933), 226-7. These 
notes will be referred to as I and II respectively. 

t Cf. Bolza, loc. cit., p. 245. The functions g and F, are related by the 
equation g = F"**F,. 
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OF (x, p) 
opt 


where | 
For a fixed value of x 
E(x,p,é)>0  (€ £ pp', 
if and only if the loci in the &-space given by 
F(x,g)=r (r>0) 
are ovaloid with respect to the straight lines.* 

First observe that these loci are homeomorphic to an (n—1)-sphere. 
For it follows from the above conditions on F that each ray issuing 
from the origin meets each of them just once. Also, because of the 
homogeneity in p and € of the functions # and F, it is sufficient to 
prove this for a particular value of r, and for values of p satisfying 
(2.2) with € = p. 

The locus (2.2) is ovaloid if and only if the tangent hyperplane 
at each point lies entirely outside the surface. That is to say, if and 
only if , ‘ 

J F(«,é)—r > 0 (2.3) 
for every value of € satisfying 

Fé‘ —p') = 0 (€' #~ p’), 


p being an arbitrary point on (2.2). But 


F,(é'—p!) = Fyéi—r 
= F(x, é)—r—E (a, p, &). 


Therefore the conditions (2.3) and (2.4) are equivalent to (2.1), and 


the theorem is established. 
Since g + 0 the quadratic form in &, 
Gis; p)E'E?, 
is either positive definite or takes on negative values. In the latter 
case it will be negative for at least one value of & satisfying the 
relation ej 
‘ Gij(v, p)p's! = 9. 
7 . , i j ~s 
For Jij(@, p)p'p’ > 9, 
and if a non-degenerate quadratic form, 
a;; ee’, 
is not positive definite, the polar plane of a given point p, 
a, pe! = 0, 


* In the two-dimensional case this has been proved by C. Carathéodory, 
Math. Annalen, 62 (1906), 449, § 2. 
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contains at least one point for which a;;é‘¢' differs in sign from 
a;;p'p', provided the latter is not zero.* 

Therefore another necessary and sufficient condition for the locus 
(2.2), which is equally well defined by 


Ji;(*, HPF = +, 


to be ovaloid is thatt 
Gi(x, p)E'E! > 0 (p Fx 0, é x 0). (2.5) 
Hence (2.5) is equivalent to the condition (2.1). Moreover, this in 
turn is equivalent to the condition 
Q(x, p,€) = Fi EE >O (EF pp’), (2.6) 
0? F (x, p) 
ep'ep) ~ 
For 915 = PRyyst Fy Fs 
and Fi = F-'9;, p*. 


wT 


where F ips — 


Hence Prins = FAG; 9ia—Iir Gd PP’; 
and it follows from the algebraic theory of quadratic forms that the 
condition (2.6) is equivalent to (2.5), and therefore tot (2.1). 

The condition (2.5) is a consequence of (1.1). For if the latter are 
satisfied, the locus (2.2) is closed and the quadratic form (2.5) is 
positive definite near at least one point p (e.g. a point on (2.2) at 
which é' has a maximum). Therefore (2.5) is satisfied since the 
signature can only vary as (p',...,p”) pass through a set of values 
for which g = 0. 


3. Just as a positive Riemannian metric determines a Euclidean 
geometry in the tangent spaces, so the more general metric 


ds = F(x,dx) 
determines a metric geometry with the distance given by 


in the tangent space at each point g. According to the last section 


* This can be seen at once by transforming to coordinates, 7, in which the 
components of the vector p are (1, 0,..., 0) and in which 
ay FE) = +yitot-..t- 
+ See I, § 4. 
t For n = 3 the conditions (2.1) and (2.6) were shown to be equivalent by 
W. Behagel, Math. Annalen, 73 (1913), 596-9. Another proof, and an exten- 
sion to n > 3, appears in the paper by Bliss referred to above. 
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a differential metric space* which satisfies the condition (1.1) is one 
which determines a metric with ovaloid spheres in each tangent 
space, a sphere being the locus of points equidistant from a given 
point. Thus in each tangent space there is a geometry of the kind 
first studied by Minkowski.+ As in Riemannian geometry, normal 
coordinates at a point qg are a means of setting up a relation of 
closest contact between the underlying space and the metric tangent 
space at g. We shall use this relation to carry over our geometrical 
interpretation of regularity from the tangent space to the underlying 
space. 
The geodesics are given by equations of the formt 

a 

Ox (2), (3.1) 

ds? | ds 
where the functions H are homogeneous of the second degree in dx, 


and the parameter s measures the are length. That is to say, 


| dz 
Fix, l. 
(2 7 
The solutions to (3.1) are given by 


x' = ¢'(q, p,8), 


( #'@-P. 0) q' 
where dd 
L (ar)... =" 


ds 


and because of the homogeneity of H they may be written as 


at = f*‘(¢,y), 

with y* = p's. 

If s is to be the are length, p must be a unit vector. That is to say, 
F(q, p) ie (3.5) 


The equations (3.3) satisfy the conditions 


6) 
oy’ 


* For an account of these geometrical ideas, introduced by P. Finsler, see 
L. Berwald, Atti del Congresso Internazionale dei Matematici (Bologna, 1928) 
263-70. 

+ ‘Theorie der konvexen Kérper, insbesondere Begriindung ihres Ober- 
flichenbegriffs’: Gesammelte Abhandlungen (Leipzig, 1911), vol. ii. 

t See M. Mason and G. A. Bliss, Trans. American Math. Soc. 9 (1908), 443. 
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and determine a transformation to normal coordinates y in which 
the geodesics through g are given by (3.4). The functions H and their 
derivatives with respect to 2 and dx are continuous for all values of 
x in the region of definition and for all non-zero vectors dx. There- 
fore the transformation to the coordinates y is defined for all values 
of x near q. 

The locus given by (3.4) with s = r and p given by (3.5) will be 
called the sphere of radius r and centre qg. It is given implicitly by 


the equation Fqy) =r (3.7) 
Let « be a vector tangent to (3.7) at the point (p's,...,p"r). That 


is to say, = 
} Fiat = 


Further, let y' = y'(t) 


be the equation of the geodesic through (p'r,...,p"7) in the direction 


x, t being the linear function of the are length determined by the 


y'(0) = p'r 
(8) -» 

dt } 
Then, as in I and II, we have 


F{q,y()}—r = $0 [Big alod +Fe H"(pr, «) ep, + OP) 


conditions 


= 1 Fy ato) + Fi | + O0(é) 
. 
1 . 

= 1{-Qcp. a+ Fy H'| +0). 


If r is sufficiently small the coefficient of ¢ has the same sign as 
Q(q.p,x). Therefore there is a positive function 8(r), such that 

Fig. y()}—r 
has the same sign as Q(q, p, x), provided 

—a(r) << t < &(r). 
Therefore the argument used in I shows that small spheres are 
ovaloid if the metric is regular. 
If the metric is not regular there is a pair of vectors, p and &, 

associated with a point q, such that 

Q(9,p.€) < 9. 
Since Fiipp' = 0, 


it follows that O(q, p.E+pp) = Q(¢. p. €) 
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for any value of p. Therefore the vector 
at = £'— pi PF €//F(q, p) 
satisfies the conditions ( Q(q.p.x) <0 
| Fi x" 0. 

Therefore regularity is necessary as well as sufficient for convexity, 
and we have the theorem: 

In a differential metric space, sufficiently small spheres are ovaloid 
with respect to the geodesics if and only if the metric is regular. 











INTEGRALS FOR THE PRODUCT OF TWO BESSEL 
FUNCTIONS (II) 


By A. L. DIXON (Ozford) and W. L. FERRAR (Ozford) 


[Received 3 October 1933] 


THE present paper contains results complementary to the 
formulae obtained in our previous paper* on the same topic. Parti- 
cular examples of our new results occur in an earlier investigation. 
We ry "a integral expressions for the products K,,(X)K,(x) 
K(X) L(a 


X)K,(x); the most noteworthy are those given in 
5 (5 ) and ” 9 which reduce when X = x to 


K,(x)L,(2) = [ Jp4)(2asinh the-H" dt. 
1 


1.2. Notation 


Throughout the work the numbers X, 2 are to be considered as 
positive unless a statement to the contrary is made. If y is any 
positive number, arguments are fixed by the conventions 
argy=0, arg(yi)= 37, = arg(—yt) = —}. 
In order to set out the work as shortly as possible we standardize 
certain abbreviations (taking the positive square root in all cases) 





Y X + ae?4i \ Mu +») 
= X2 2 2 3 s 2 wi —_ —2v0i a 
A = .(X?+2?+2X-2 cos 26), e € (x=) 


er . X + xe-2"\ Ku +”) 
A = ,/(X?+a?+ 2X2 cosh 2u), @ =z com 


X +e?" ; 
a = }log(X/x); 
when X >a and1<e"“< X/x 
] on |X —xe-2™\ He +) 
Ay = .(X?+2?— 2X2 cosh 2u), ew = eo : 


X —ae2" f 
when X > 2 and e®" > X/x 


A, = 4/(2X2 cosh 2u— X2—2? —_—_—— 
2 v( ), wve2"— X 


Quart. J. of Math. (Oxford), 4 (1933), 193-208. References to this paper 
are denoted by I. 
+ Ibid. 1 (1930), 122-45; at 143-5. 


y —are-2u\ Mn+) 
et: — ens ve ) . 
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2. We begin by noting that the formula I (6.35), wherein 
= ,/(X?+2?—2X-2 cos 9), Q = .(X?+2a?+2Xzcosh wu), 
may be written as 


27 sin(u-+-v)a J, (XJ, (2) [X >a, R(p—v) < 3] 


»(a)} dp — 


{sin pr J), (w)+sin ve J_,, 


(Q)} du. 


. = =I. re" d(+v) 
2 sin vr fess {sin war sin yr 
—2sin va e- ( {sin yar J), ,,(Q)-+-sin va J_, 


X+axe-" 
0 

In this form, when Y,(X) replaces J,(X) on the left-hand side, 
the only changes on the right-hand side are that the functions Y,,,,, 
Y_,-, replace J,,,, J_,-,: the statement is proved by calculations 
similar to those whereby I (6.35) was established. 

3.1. We next set out the key transformations for our remaining 
formulae. We assume that | R(u)| < } and | R(v)| < }, and we leave 
aside the question of analytic continuation to other values of » and v. 
Further, we suppose that X > 2 > 0 unless the contrary is stated. 


Consider 


X —wxe2U B+ 


- yop —2U \ +) 
aad = )" K,,.f\(XeU —ae-")(xe’ —Xe-V)} dU, (1) 


0 
where the path of integration is the real axis indented upwards at 
the point U x, and the conventions of I (4.41) apply to the 
multiple-valued functions. 

'7 and use the substitutions 


When we increase the argument by 3} 
of I, § 5, we see that (1) is equal to* 
ba 0 
erm’ | ebik 4 | evmi | ewk ; 
ve} | ew K, . (td) dé +-er7' eK, ,,(tA) du. 
0 0 
In later stages it is easier to deal with multiples of exp($77) if we 
write (1) as 
x = 2 


| eK, (tA) du +eitH+) eK 


pty pet 


(Ag) du. (3) 


0 x 


* The result I (5.22) is, in fact, obtained by adding the present result with 
p, v to the corresponding result with —p, —v. 
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The equivalence of (1) and (2) may then be set out as 


[er Ky siMy) du Helmer) [ emeK, (Ag) du 
in . ea) 
= ier | eK, ,,(id) do +e" | eK, ,,(iA) du. (4) 
Q F 


The equation (4) has been proved on the supposition that 
X >a > 0. We may think of it as a formula which involves the 
two variables X and 2/X (= €, say). Having written it in terms of 
these two variables and keeping € fixed, we find that, by analytic 
continuation, the equation is true for all X in the half-plane R(X) > 0 
and that, by continuity, it remains true when arg X = —4z (though 
not when argX = $m), provided always that R(u+v) <1 and 
R(v—p) < 4. The latter conditions are satisfied in view of our 
hypotheses | R(u)| < 4, |R(v)| < F. 

We may, then, write X exp(— 4?) for X and xexp(— 477?) for x in 
(4) and so obtain the formula, valid when X > x > 0, 

[ eK, 4(Aq) du + ebminsw) [ eK, 44(—iAg) du 


0 
i a) 


har 
= ie [ eK, (0) dd +e" [ eK, .,(A) du. (5) 
0 0 
For brevity we write this formula in the notation 


O(y,v,i) = ie'f(y, v, i) +e" F (u,v). (6) 
This is the key formula to much of our later work. 


3.2. We may, of course, use the original form (1) and write 


ao 


. X —xe-2"\ i+») le 
be nae | em X —xe2" Ky avy (X —xe®)(X —xe-*")} du, 
0 


where the arguments of the factors of the integrand are zero for 
» = 0, and at u = «the path of integration has an indentation above 
the real axis; or in a form more amenable to certain applications, 


- —a/t\ie+) ; 
D(y,v,i) = t [ e(5-3) = K,,.,{4(X—2/t)(X—at)} dt, 


where arguments are zero for ¢ = 1 and the path of integration has 
an indentation above the real axis at t = X/zx. 
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It follows that ®(u,v, —7) is defined by an integral whose path has 
an indentation below the real axis at / = X/z. 


4.1. Identities 
The identity (6) of 3.1, namely, 

D(u,v, i) = ie’**f(p, v, i) +e" F (u,v), (1) 
shows that any product of Bessel functions which can be expressed 
in terms of the ® type of integral can also be expressed in terms of 
the f and F types of integral. This leads to a duplication of formulae 
which, when encountered in isolated examples, is rather puzzling. 

There is a further identity, involving only ®-integrals. This follows 
from (1) and the fact that 
f(u.v, —1) = f(—p, —v, 1), (2) 
which is evident since K,,,,(v) = K_,_,(). We have, from (1) and 
its conjugate complex, 
e7'D(u, v,i)—e” "O(n, v, —2) = if f(u,v, 2) +f(u, v, —d)}, 
and this, by (2), is ie htA—p, 
But the latter is unaltered if we write —yu, —v for p, v, and so 
e™D(u, v, 1) —e”""O(u, v, —7) 
= e'"M(— pn, —v,1)—e-""'O(— p, —v, —7). (3) 
This identity is of service in reducing later results to simpler forms. 


4.2. Before using the identity (3) we show that it is, essentially, 


very simple in character. It can, with some labour, be proved as 


a direct consequence of Cauchy’s theorem. 
For large values of |t} and for the arguments of X—at, X—z2/t 
which we shall consider, 
_1[X—a/t\ie |, Bw " 
p aa] Ku sty(X —at)(X—x t)} 
is at most O(|t|/!”-4#-#) = o(|t|-1), and so, by Cauchy’s theorem, 
(0+) , 
ry ([X—2/h\iem |, . ae 
0 | t (3) Ki sty(X —at)(X —x/t)} dt, (1) 
where the contour has indentations at the points t = X/x, t = x/X, 
and argt = 0, arg(tX —x) = 0, arg(X—at) = —7 


at the starting-point. 
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This is the identity (3) of 4.1. On taking the parts separately, and 
using the form given in 3.2, we find 


from 0 to 1 contributes —2®(y, v, 7), 
from 1 to0 contributes —20(—y, —v, —i), 
and, on the return along the lower side of the real axis, 
from 0 to 1 contributes 2¢?7"'®(—y, —v, i), 
from 1 toc contributes 2¢e?7"'@(y, v, —i). 
5. Expressions for K,(X)K,(x), K,(X)/,(x) as integrals of the 
® type; X > x 
The formula I (4.42) may be written as 


a 


K,,(iX)K,(—ix) = | (e+e) K,, (Ay) du + 
0 


elm { eK, (Mg) du ern | e-sK,, (Me) du. 


a a 


In this we may put —iX, —ix for X, x and so obtain 
K, Ea (xe~**) ax Pigs t)+@(—p, —v, 1%), 
ie. eK, (X)K,(x) +i K(X = D(y, v,7)+0(—p, —v,i). (1) 
From (1) and its ai: ge we have 
2 cos vm K,,(X)K,(x) 
(u,v, 1) +O(—p, —v, —t)+O(p, v, —1)+0(—p, —v, —1), (2) 
or, on using the identity (3) of 4.1, 
K,(X)K,(«) = e'*{®(—p, —v, i) + Ou, », —i)} 
= e"7{@(p, v, t)+0(—p, —v, —1)}. 
We also have, from (1) and its conjugate complex, 
27i cosvm K,(X)1,(x) 
= eY*{O(y,v, i) +O(—p, —v, )—e™{O(y, v, —i) + O(—p, —v, —i)}, 
or, on using the identity (3) of 4.1, 
mi K,(X)1,(x) = O(—p, —v, 1)—O(—p, —v, —1). (4) 


The equation (4) reduces, almost at once, to the curious formula 


ca 
) = fem Iera a) du 
a 
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or, in full, 


K,(X)L,(2) 


o 4) 
. F X — xe-2"\ —Me+) ; - : " 
= | e “ Jsyly(2Xa cosh 2u—a*—X*)} du. (5 


} xe?"— X 
x 


The simplest form of (3) is then seen to be 
2i sin v7 K,(X)K,(x) 
= O(—p, —v, —1)—®(—p, —v,1)+- O(y, v, 1)—O(p, v, —2), 


or 


2 sin vr K,(X)K,(x) 


0 
* oy, _X —xe-2u\ e+) ; ae 
= 7 | en J_, Ay (2X cosh 2u—x?—X?)} du — 


we2"— X 


oo) 
P y —ge-2u —h(p+v) 
o—2vu f 9V ranch 97, — v2 72 > 
—T | "aie ( = c) Jiri (2Xx cosh 2u—x*—X?)} du. (6) 


6. The product 1/,(X)K,(x);X > «x 

The task of finding a compact formula for /,(X)K,(x) is by no 
means a simple one. We begin, as we did for K,,(X)J,(~), by quoting 
I (4.42) in the form 


x 


K,,(iX)K,( —ix) = [ (e+ e-) K,, ,,(vAy) du + 
0 


+ eizi(uty) | evr K., y(Ag) du +e-t7H+) | e-vs Ky, +: 


A.) du. (1) 
x x 
In this we may put 7X, ix for X, x and so get 


x 


K,,(Xe™)K,(x) = | (e+e) K 


u+v(Aye™) du + 
0 
a a) 
+ ebti(u v) els K,, (iAg) du 4 e— 7i(u+r) ens K 


x x 


We now note that, from the definition of ® in 3.1 (6), 


f enek 


utvltAg) du. (2) 


(tA,) du = b+ O(n, v, —t)— | eK, (Ay) du}, 


a 0 


perv 
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and, on substituting for the «, co integrals, we write (2) as 


eu! K (X)K,(x)—mi 1,(X)K,(x) +2 | cosh wy J, ,,(A,) du 


0 


a 
= eter f (ee) Kgl My) de + 
0 


eon iy, v, —i)— [ ets Ky +(Ay) du} + 


0 
a 


fewer 6p, —vy, —t)— eK, 4 (Ay) du} 


0 


—2¢sin(u+v)z [ e"1 Ky +(Ay) du + 
0 
+eHmri@(n, v, —i)+eH"'O(—p, —v, —i). 

Now transfer the K,(X)K,(x) term to the other side of the equation 
and use the sconiiadin [§ 5 (3)] 

K,(X)K,(x) = e*™{(y, v, i) +O(—p, —v, —i)}. 
This gives us ‘ 

—at 1, (X)K,(x)+-2m0 | cosh wy L,4,(Ay) du 


0 
a 


—2isin(u+-v)z | e" Ky (Ay) du + 
. +e! O(n, v, —i)—e-“H""'@(p, v, 2). 
When we insert the values of the functions ®, this becomes 
eluty)at Fen Ky 4,(tAg) du —eKu +r fen K 
or, on using 1 (6.13) and (6.14), : 


—ai [ ered. (Ag) du. 


(—7A,) du, 


pty 


We have thus proved that 
1,(X)K,(a) = 2 | cosh w, J, ,,(Ay) du + [ ed, .4(Ag) du. (3) 
0 a 
7. Expressions for K,(X)K,(x), K,(X)J,(x) as integrals of the 
fand F types; X > ax 
When we make the substitutions 


d(u, v,t) = ett f(u, v,)+F(p, v)}, 
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and so on, in formula (2) of § 5 we find that our result reduces to 
the integral which formed the starting-point of our investigations in 
I. Thus our I (2.5), i.e. 


K,(X)K,(x) = | e”K,,,(A) du, 
is the f, F form corresponding to the formulae of § 5. 
The formula (4) of § 5 becomes 
m1 K(X) L(x) = e7i{a f(—p, —v, 1) + F(—p, —v)}— 
—er7{__¢ f(—p, —v, —t) ++ F(—p, —v)}, 
and this, with a little simplification, gives the result 
2a K,,(X)L,(2) 


; 4 phi \ Mu+r) 
-- | ¢ “4 : = ai} , Kivi (X?+2?— 2X2 cos $)} dp — 
A —xe-?" 


= 


wo 
7 


X +2xe- u\ d(+v) ’ 
—2sin vz | er( Ys :) Ki svty(X?+2?+ 2X2 cosh u)} du. 
; X+2xe 


This formula is the analogue of those given in § 2. It can, in fact, 
be checked by using the formulae of § 2 to derive integrals for 


2a K,,(—iX)J,(x)e byt 








A REDUCIBLE CASE OF THE FOURTH TYPE OF 
APPELL’S HYPERGEOMETRIC FUNCTIONS 
OF TWO VARIABLES 


By W. N. BAILEY (Manchester) 


[Received 19 September 1933] 


1. As long ago as 1882 Appell* proved the formula 
. X—Y 
aa Bs; =) 


(1.1) 


where a; B, B’; yi X, Y) >2% (x) moll) n(B) n xm Y*, 
m ned n= Y)m+n m! n. 


thus showing that the function F, can be expressed in terms of an 
ordinary hypergeometric function when y = B-+-f’. 

No formula of such generality was obtained for the fourth type of 
function, namely, 


F(a, B; y,y'; X, Y) = S > & (24)m+n(B)m-+n xmyn, 


A, A (Yl ¥')n min! 
although Appell showedf that the functions 
Fila, B;y,a+B—y+1;2%,(1—z)?] and {F(«,B;y;2)}? 
satisfy the same differential equation of the third order. 
Some years ago Watsont{ gave the formula 
PoP BLO —a— BF Bs 70+ B—y+ 1322, (1—2)(1— 2) 
+{(1—2)(1— Z)}7-* PD (y—a) (yy) (a+ B—y) X 
x Kly—«, y—B; y, y—a—B+ 1; 2Z, (1—z)(1—4)] 
= D(@)P(B)P(y—a)P(y—B) F(a, Bs 32) F(a Biy;Z), (1.2) 
and he stated that this appeared to be the best theorem concerning 
the expression of functions of the fourth type in terms of hyper- 


geometric functions.§ 
In this paper I give a simpler theorem of this nature, in which only 


* Appell (1), 175. 

+ Appell (2), 418-21. Appell stated the result explicitly only in the case 
when y’ = a+B—y+1. t Watson (5). 

{ Watson also gave a relation connecting four series F, with two products 
of hypergeometric functions, but this can easily be derived from the formula 
quoted above. 

3695-4 x 








306 W. N. BAILEY 
one series F, appears. Watson’s formula is really a combination of 
two formulae of the type given in § 2. 
2. The theorem is that 
Fa, B; y,«+B—y+1;2(1—Z), Z(1—z)] 

= F(a, B; y;z) F(a, B;«+B—y+1; Z), (2.1) 
this formula being valid inside simply connected regions surrounding 
z = 0, Z = 0 for which 

|2(1—Z)|#+ |Z(1—z)|# < 1. 

The theorem is true, for example, if0 << Z< land0 <z< 1—Z. 

This result shows that the function F, can be expressed as a product 
of ordinary hypergeometric functions when y+y' = «+B+1. 

It will be noticed that, when Z = 1—z, the function on the left of 
(2.1) is that considered by Appell, and the function on the right 
satisfies the same linear differential equation as that satisfied by 
{ F(x, B; y;z)}*. 

The theorem can be proved by an argument similar to that adopted 
by Watson, but in this proof it appears to be necessary to use the 
formula* — 

C D(ay-+s)I(a+8)I(a9-+8)F'(1—B,—8)'(—8) ds 
I'(B.+8) 


I(a,)10( ¥5) (a3) (1 — B+ x, ) 00 —B,+ X)*(1—B, +3) 
TP'(B.—24)1"(B,—ag)1"(B,—avg) 


where By +B. = %+a.+0,+1, 


instead of the more familiar Barnes’s lemma used by Watson. Con- 
sequently, I give the more elementary argument which now follows: 


F[x, 8; y,«+B—y+1;2(1—Z), Z(1—z)] 


(X) m+n(B)m-+n {2(1— Z)}"{Z(1—z)}" 


OTE ye 


(Y)m(a+B—y+1), m!n! 


(x) \(B)~ f2( on Z)}"{Z(1—z)}N-™ 


(Y) nf x+B—y+ 1)y-m m!(N—m)! 


( x) (B)y ZN(1—2)F|”— a—B—N, —N. (z =| 


Ni(a+B—y+1)y y 1-2 


N=0 


* Barnes (4). As usual with integrals of this type, the contour is curved, 
if necessary, to separate the increasing and decreasing sequences of poles. 
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_~ 5 ww gx _ay-0-8 
- 2. 
x F[y—a—B—N; —N, y+N; y;2/Z, z] 


ax) (8B) Z* 
(1—z)r- a> > > 


=0 m=0 n= 
,y—a—B—F),..l— N) n(y+N MG zy zn 
In! Z : 
(Y) m+n m:n. 
Here we have used (1.1). Now put VN = p+m, n = q—m, and we 
obtain 
(1—2)7-2-B x 
B) »n+m(¥— a— -B)q- ~p—m' (y+p+m),- m ZP2t 


«> > 2! ca Eleent P\(y)q m! '(q—m)! 


ap SN (p(B —a—f) (Y+P)q 2 
— er Pig'(y), er 


* ed B+p, —q | 
viii ytp, l+a+p—y—q+p 


The series ,F, can be summed by Saalschiitz’s theorem, and we 


obtain 


pZ” — (y—)gy—B)_2* 
—z)v-«-B qi lq 
a > > 7 p\( me 7, Vy)q 


p=0 g=0 


s (1—2y-+-##| B; “| ie 


atBp—y+l Y 


=P apt 
‘a, a+B—y+1 
In this proof we assume that |z| and |Z| are sufficiently small,* and 
then the complete result follows by analytic continuation. 
3. Now in (2.1) change z, Z into 1— Z, 1—z, and we find that 
Fl x, B; y,«+B—y+1;2(1—Z), Z(1—z)] 
= F(a,B;y;1—Z)F(a,B;a+B—y+1;1—z) (3.1) 
inside simply connected regions surrounding z = 1, Z = 1, for which 
\Z(1—2)|*+ |2(1—Z) |! <1. 
In particular it is true if 0 < Z< land1—-Z<z< 1. 
This result, together with (2.1), gives the complete expression of 


* It is sufficient to assume that |z] and |Z| are each less than }(v2—1). 
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the function F, in terms of ordinary hypergeometric functions when 
yty’ = a+B+1. 

4. It is perhaps worth while to justify the assertion, made at the 
end of § 1, that Watson’s formula can be derived from (2.1). For, 
using the theorem twice, we see that the left-hand side of (1.2) is 
D(a) E(B) P(y)P(y—a—B) F(x, B; y3 2) F(a, B3 a+B—y+1;1—Z)+ 

+{(1—z)(1—Z)}7-* FP (y—a) P(y—B)P'(a+B—y) Ty) x 

< F(y—a, y—B; y32z)F(y—a, y—B; y—a—B+1; 1—Z) 

= F(a, B; y; z)[T( x) P(B)P(y)P(y—a—B) F(a, 8B; «a+B—y+1;1—Z)+ 
+T(y—a)P(y—B)P'(a+B—y)P(y) x 

x (1—Z)r- BF (y—a, y—B; y—a—B+1; 1—Z)] 
> F(x, B; y;2)P(a)P (8) (y—a) 0 (y—B) F(a, B; y; Z). 

This argument assumes that z, Z lie in certain regions surrounding 

z= 0, Z = 1, but the symmetry of the result shows that it is also 
true in the regions surrounding z = 1, Z = 0. 

5. Finally I mention that the present investigation of (2.1) was 
suggested by the formula* 


F(a, B; y3z)F (a, B; a+B—y+1;z) 


= nfrF Hoth), Hot Br atet—a)) 
y,a+B—y+1, a+f8 

The reader will readily see the connexion between this result and 

the special case of (2.1) in which Z = z. It can, in fact, be derived 


from (2.1) by a use of Vandermonde’s theorem. 
* Bailey (3), § 6. 
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ON A PROBLEM IN THE ADDITIVE 
THEORY OF NUMBERS (VI) 


By ©. J. A. EVELYN (Wotton) and E. H. LINFOOT (Bristol) 


[Received 6 September 1933] 
1. Introductory 
In previous papers we discussed the representation of a large number 
n as the sum of s Nth power-free numbers or M-numbers, that is, of 
s numbers not divisible by an Nth power greater than 1. Here N is 
fixed once and for all, repetitions are allowed, and order is regarded 
as relevant. Using an elementary argument, we were able to showt 
that for N > 2, s > 2 the number of representations 
s-1 = 
v,(n) = = De a S(N,s8;n) [1+0(n™ ‘xvi *)} (1.1) 
as n > 00, where the factor S(N,s;n) lies between two positive abso- 
lute constants; the argument was later simplified by Estermann.{ 
Using analytic function theory, on the other hand, we obtained§ for 
N > 2,8 >3 the formula 
°-t 1 
vale) = yl BON) 
which is an improvement on (1.1) when s > V-+-2. 
In the present paper we extend these results to the case in which 
the M-numbers all belong to a given arithmetic progression, using an 
adaptation of Estermann’s argument, followed by an induction, to 


, 1 
1+y+ 


S(N,8;n){1+O(n 


(1.2) 


prove 
TueorEeM A,. Lets >2,N >2,1<b<a;lea=T{ p} be the 
i=1 
expression of a as a product of primes in ascending order; let (a, b) be 
Nth power-free,|| and let 


S(n) = S(a,b, N,s8;n) 


— ])s+1 —l1 
~_ 7 2S A hee 
(p r~ C... 1)8 (p} pV Og 1ly-? 
+ Evelyn and Linfoot (2) (II), (III) (see list of references). 
{t Estermann (1). 
§ Evelyn and Linfoot (2) (V). 
|| Unless (a, b) is Nth power-free, the progression (1.3) evidently cannot 
contain an M-number. 
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Then the number of representations of a large number n as the sum of 
s M-numbers, each belonging to the arithmetic progression 


b, b+a, b+ 2a,..., (1.3) 


repetitions being allowed and order being regarded as relevant, is 


1 (ay 77 LY ends ol -tepe wl 
y(n) = gayle) ia (1 ) S(n)\1+O(n N+ )| 


Dil ab) (1.4) 
for every « > 0, if n = sb (moda), while for n # sb (moda) 
y(n) = 0. 
By using the Winogradoff-Gelbcke argument of (2) (V) we are able 
to obtain the (sometimes sharper) 


THEOREM ig For s > 3, N > 2 and with the notation of Theorem A, 


rs) = e—ayila) i] (1— yx) Stop{1+-0(n +949) 


pilin b) (1.5) 
if n = sb (moda), while for n ~ sb (moda) 
y(n) = 0. 
As the proof is very similar to that for the case a = 1 we do not 
reproduce it. 
To interpret the a 


eam) TT (pes) 8 


Pilla, b) 


we observe that the number of representations of n as the sum of 


ey : ; ce : 1 n\e—1 
s positive integers of the form ar-+-b is asymptotic to @—i\e as 
, s—l)!\a 


n -> 00 through values congruent to sb (moda). To pass from this to 
the representation by M-numbers we naturally introduce a factor 


Finally there will be a certain interference, depending on the form 
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of n, through the restriction imposed by the equation 
M,+M,+...+M, =n 

on the arithmetic structure of the M’s; this corresponds to the 
factor S(n). 

Just as in the special case a = 1 of (I), we can infer at once from 
our asymptotic formula 

THEOREM B. The number of partitions of a large Nth power-free 
number M into s Nth power-free numbers of the form ar-+-b, where (a, b) 
is Nth power-free and M = sb (moda) differs asymptotically from the 


‘normal’ number 
1 M\:— 1 \s 
7 | | 1— — 
aain(s) J ( p> :) 


y a 
Pid (a,b) 


by a constant factor y(s,N,a,b) which is less than one when s is even, 
qreater than one when s is odd. 

2. We first suppose s = 2. The numbers under the summation 
signs have the following meanings: 


1. aut+ty%v = n; 2. «Nu = b (moda); 3. yXv = b (moda); 


4, (2. *) ne (v.<) = 1; 7. p(x),u(y) £0; 8. ay <nNA; 
c c 


9. (x, y)¥-*|"; 10. (a, y)*|n. 
Cc 


Si > . ° 
ince (m) = 1 if n is an M-number 
my |n ° 
= 0 otherwise, 


the number of representations 
LLY,U,0 
von) = 1 = F plx)uly) 


M,+M,=n 1.2.3. 
M,=b(moda) 


Ly uv 

= mla)uy) > 1. (2.1) 

Letc = (a,b) = Il p%, where every o; < N—1, and write (c, 2’) = 8. 
i=1 
Then c/5 is prime to x and the congruence 
xu = b (moda) 

shows that c|au, (c/8)|u. It can therefore be written 

aN «& 


FT oa b, (mod ay), 
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where a = ca,, b = ch, (a,,b,) = 1. This has solutions in wu only if 
2/8 is prime to a; it is then (writing u, = du/c) 


rN 


b, soc, (mod a,), (2.4) 


rn 


where it is immaterial which value of the socius is taken. Let u,, be 
the least positive solution in u, of this congruence; thus 1 < u,5 < a. 
Then the general solution is 

Uy = UWot70, (r > 0). 


Thus (2.2) has no solution in w unless 


aY oa 


/ a 
i.e. unless (>. i: 
Pe 


when this is satisfied its general solution is 


c ra 
= Uu i. — © 
(c, 2") 40°" (c, a ) 


A similar result holds for v. fig therefore becomes 


where 6. stands for 


( aNr 
a 
(Cc, 2 


x 
We now introduce a new notation. Let m = [JT p‘ be any positive 
b=1 


t 


we 
integer. Then by m‘-* we shall mean the integer [] pi’-%. Since 
i=1 
v, y are quadratfrei, condition 6. can be written 
a(aN-*r+-yN—-*s) = n—ca’-*u, »—cy*—* 0, 9. 
This has a solution in 7, s only if 
tt, ; : 
.N—* 9,N—* .N—* yN-*y) 
~(2 ies U0—-Y “1,0: 
c c 
which is equivalent to the pair of conditions 
—— .N—* 9,N-«)|”% 
a|n—2b, (2*-* ¢-*)}i— 
since a/c is prime to (x\-*, yN-*) = (a, y)\-* and a4-*u, 9, y*-*v19 
are each congruent to b, (moda/c). Thus v.(n) = 0 unless these are 
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satisfied; when they hold, 


ya SS en aN—* uty o—y—*0,9)a(x, y)—* 
cy aaxN-* qyN-* 





+26; 


by Lemma 2.1 of (II), where @ denotes generally a number absolutely 
less than or equal to 1, 

_ n (2, (x, y)*-* 

~~ @ aN-#yN-* + 46, 
SINCE Uy 9, Vig < A. 


Let =F mary) 31; 


then, just as in Estermann’s original argument, 
t +€ 
\va(n)—Q,| = O(n"), 


while now 2,= 
4.5.7.8.9. 


when a|n—2b6, zero otherwise. The condition 9. (x, y)*-*|= is equi- 
c 


valent to 10. (x,y)‘|n. Thus, in virtue of 8., 


== DS ule wo) +0(n¥*), 
4.5.7.8.10. 
when a|n—2b. We next drop the condition 8. x,y < n+ in the 
first term; the error caused is O(n?(*+) as with Estermann. Thus 


for a|n—2b 
ru 


_% (x, y)>~ + 
v9(n) = ~ ‘ Mea Y) ye a" O(n > 
-5.7.10. 
while v,(n) = 0 if aJfn—2b. 
The coefficient is easily evaluated; by a calculation similar to that 


of (2) (IIL) we readily verify that for s > 1 


H(%)  p(2,) r 
_N-e os (Xqy0005 x,)" " 
8 


a) TT (+o sap) IT (+ epee 


N a N a 
pi In, rif D; In, it, 


) S(a,b,s,.N;n). 


Taking.s = 2, we obtain Theorem A, for the case s = 
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Next let s > 3 and suppose the theorem proved for v,_,(n). Thus, 
if a|n—(s—1)b, 


(Ly,.009 Ls-1) ¥ | n 


for every « > 0, while if a|n—(s—1)b, v,_,(n) = 0. Whence 
n 

m= S 1= YF ya(n—M,) 

M M,=n 


ites. P M,=1 
M;,=b(mod a) M,=b(mod a) 


n s—2 Tiyrees Ls 

(n—M,) . 

7 D)\' 48-2 

8—<a) a” " 
M,;=1 ( ) 

M=b(mod a) 

(Xj 50009 Ls-1)¥ | n— Ms 


n s—3+ 2 +e 

+O 3} (n—M,) “N+i", 

Mi=1 
M,=b(mod a) 

provided a|n—M,—(s—1)b, i.e. a|n—sb, while v,(n) = 0 if a|n—sb. 
In the former case the rest of the argument is similar to that on 
pp. 641-4 of (2) (III); we therefore merely state its result, namely, 
that for a|n—sb 


] 
a (s— a 


(Lry0005 Lg) jn 


while for a/n—sb, v,(n) = 0. In virtue of (2.5), this proves Theorem 
A, by induction.t 


+ On p. 641, 1. 16, of (2) (IIT) an error term O(n aaa *) is omitted. 
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ON A MEAN FREE PATH FORMULA 
By E. A. MILNE (Ozford) 
[Received 17 May 1933] 

1. Tue following note arose out of the recent brilliant work of 
Blackett and Occhialini* disclosing the probability of the existence 
of a positive electron. Following the ideas of Dirac, Blackett and 
Occhialini evaluated the probability that a positive electron will dis- 
appear by reacting with a negative electron to form two photons 
during a length of path in which its energy, if it had not been anni- 
hilated, would have fallen from £, to #,. A formula is obtained 
below which relates this probability to the mean free path for 
annihilation. The analysis differs from the usual analysis in the 
kinetic theory of gases owing to the fact that the future mean free 
path of a particle is a function of its present energy, and so varies 
along the path; in the kinetic theory of gases the future mean free 
path for collision (for a particle of constant velocity) is independent 
of the portion of path already traversed. 

2. Let f(£)e be the probability that a particle is annihilated whilst 
its energy, if it had not been annihilated, would have fallen from 
E to E—e (e small and positive). Let dr be the corresponding length 
of path for non-annihilation, and write 

dr 


an —9(#) (g(£) > 0). (1) 


Let ¢(£,, £) be the probability that a particle avoids annihilation 
whilst its energy falls from 2, to E, r being the corresponding length 
of path described. 

Then if N particles start with energy L,, N¢(#,, £) survive to 
energy EF, and N¢(#,, Z—e«) to energy E—e. The loss of particles 
between F and H—e is accordingly 

,2(Ey, E) 
N———«, 
ok 
and this must be equal to 
N¢(Ey, B)f(B)e. 
6¢(E,, E) 
= ¢(F,, E)f(£), 
* Proc. Roy. Soc. A, 139 (1933), 699. 


Hence 
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Ey 
or log $(EZ,,E) = — | f(E)dE (BE < #)). (2) 
E 
The particles annihilated in the interval (7, H—e«) have a range r. 
Hence the mean range, or mean free path for annihilation, of particles 
starting with energy H,, say A(F,), is given by 


E=E, 
‘ E : E : 
BE, = [ Ae), ap. 
. ok 
E=0 
Integrating by parts and noting that 

at EH = E,, r= 0, 

at HE = 0, ¢(£,, Z) = 0 


(assuming that all the particles are annihilated by the time EZ = 0), 


we have E, 
. , dr 
ME,) = — | $(B, 2) 7 a 
0 
Ey 
= | $(H,, E)g(E) dE, (3) 
0 
by (1). Hence, by (2), 
ae 
rp» — | f(E) dE . 
NH) = [et “4(B) a. (4) 


0 
We now differentiate with respect to H,. The details are most con- 
I 1 
veniently carried out by writing 


E, 
f(E) = FB), — | f(b) dE = F(B,)—F(2). 
E 
We find . 
7 rs ~ {sue ae 
dM Ey) = 1) | ez g(E) dE 
dE, 
= 9( 7) —f(B,)ME,), (5) 
by (4); (5) is the essential formula. Omitting now the suffix 1, we have 
{(B) = g(E) 1 ddA(B) 


NE) XE) dE 
_ 1d _1 dX) 
N£)dE XE) dE 
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Inserting in (2) we find 
=0 
E. 1) 
log (E,, E) = log 
r=r(E;,E) 
E =E,,r=r(E,,E,) 


d NE, 
log $(#,, E,) = — sant BNE 


E=E,,r=0 
In Blackett and Occhialini’s notation, 
¢(H,, H,) = 1—®( E,, E,), 
so that (8) may be written 


log[ 1—®(£,, E,)| ~~ fi + logs 


aE ) 
(Ey) 
This is the desired formula. 


3. Introduce now a new variable R defined by 
(Fy, E) — R(E,, E)—X(£), 
so that* dr = dR—X'(E) dE. 
Then (9) reduces to 
E; 
dR 
log[ 1—®(£,, E,)| = — a 
og[ O( 1» 2)] X( Ey’ 
Ey 
where at E ie E,, R ad \(E,), 
at EK = E,, R = X£,)+7(£,, £,). 

The physical meaning of R(E,, Z) is that it is the total path which 
a particle would describe if it (a) escaped annihilation whilst de- 
creasing in energy from E, to EF, (6) then pursued a path equal to 
the mean free path for annihilation from energy ZH. It may be noted 
from (10’) that dR is not equal to the element of path dr unless 
\'(£) = 0, i.e. unless the future mean free path for annihilation is 
independent of the energy. 

4. In the elementary treatment of mean free paths in gas-theory, 
the chance of a collision within a further distance dr is simply dr/A, 
where A is independent of the distance already traversed. This agrees 
with (5) on taking A(#) independent of Z, since the chance in question 
is —f(E£)dE, ie. —g(E)dE/A, i.e. dr/A. The factor A(£,)/A(Z) in (9) 
then, of course, reduces to unity. 

* dR depends only on E and dE. 
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5. Writing (8) in the form 
7 i dr 

48, 8) MEd h as 
we readily verify that the right-hand side of (12) is necessarily less 
than unity, for #, < #,. Formula (12) is fully analogous to the 
formula for the density p at height r in a non-isothermal atmosphere, 
namely, 


Po P 
\(£) playing the part of temperature 7’. 


7 mg fdr 
Pp 1 k j T 


6. I am indebted to Dr. Blackett and Professor Dirac for pointing 
out that in the formula used by Blackett and Occhialini,* namely, 
E, 
- hdr 


log[ 1—®(£,, E,)| = — \(B)’ (13) 


ss = 


E, 
A(£) is not the actual mean free path of a particle of initial energy 
E but the mean free path for annihilation assuming the particle 
continued with constant energy /, as would be the case, for example, 
if a suitable external field were applied; (13) is equivalent to (2), 
since —f(#) dE is here dr/A. It follows by comparison of (13) and 
(9), or by direct calculation of A(#), that 
\(E£) ogg Ge 
1—(B) sea. dr’ 
which is the relation between A(#) and A(#). Naturally always 
A(E) > X(£). In Blackett and Occhialini’s paper r and A(£) are 
tabulated numerically as functions of 2, and ® is computed from (13). 


(14) 





* Loc. cit. 715. A misprint of sign has been corrected. 














THE UNIQUE FACTORIZATION OF A POSITIVE ” 
INTEGER 


By F. A. LINDEMANN (Ozford) 
[Received 3 August 1933] 


THERE are several demonstrations of the theorem that it is not 
possible to split up a number into prime factors in various different 
ways. The following very brief and simple proof appears not pre- 
viously to have been published. \ 

The argument rests upon the fact that it may be shown that if 
there were any number which could be resolved into primes in more 
than one way, one could always find a smaller number of which the 
same was true. By repeating this process, we obtain a decreasing 
sequence of positive integers each of which can be resolved into 
primes in more than one way. But such a sequence must eventually 
contain a prime or unity, for either of which decomposition is unique. 
We have then a contradiction. 

Suppose there is a number ” which has two different compositions 


N = Pi P2Ps--- = 19293 I--: « 
We wish to show that in this case there must be a smaller number 
n’ with the same property. 

We may suppose that none of the factors p,, po, ps is identical with 
any of the factors q,, %2, 73,--., for obviously if any q were equal to 
any p, one could cancel on both sides and thus obtain a number n’ 
smaller than the original number n having two different decomposi- 
tions. Further, we may suppose that p, is the smallest of the p’s and 
q, the least of the qg’s. Then one of these must be < vn and the other 
< vn, so that p,q, <n. 

Hence Pr P2Ps— Ph = 192%39e-—Pih 
is a positive number and 

Py(P2P3---—G1) = 91(92. 93 Va--» —P1)- 

There are two possibilities: 


(a) q, divides py ps... . 
Then n’ = py, pz..., smaller by the factor p, than n, has two decom- 
positions, one of which contains g, whilst the other does not. 
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(b) q, does not divide PePs---- 


Then it does not divide p, p3...—q,. 

Hence n’ = py(Po Ps---—G1) = Py Mg Mg. 
where no « is q,, whilst since 

nm’ = (9293 I4---—P1) 
it also has a decomposition in which one prime is g,. Thus n’, which 
is smaller by p,q, than n, has two different decompositions. 

Since the same reasoning can be applied to the number n’ as to 
n, the process can be repeated as often as one likes, so that if any 
number could be split up into prime factors in two different ways, 
it would be possible to construct numbers as small as one pleased, 
which had the same property. 
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